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Abstract 

We study a nonlinear, unsteady, moving boundary, fluid-structure 
interaction (FSI) problem arising in modeling blood flow through 
elastic and viscoelastic arteries. The fluid flow, which is driven by 
the time-dependent pressure data, is governed by 2D incompressible 
Navier-Stokes equations, while the elastodynamics of the cylindrical 
wall is modeled by the ID cylindrical Koiter shell model. Two cases 
are considered: the linearly viscoelastic and the linearly elastic Koi- 
ter shell. The fluid and structure are fully coupled (2-way coupling) 
via the kinematic and dynamic lateral boundary conditions describing 
continuity of velocity (the no-slip condition), and balance of contact 
forces at the fluid-structure interface. We prove existence of weak so- 
lutions to the two FSI problems (the viscoelastic and the elastic case) 
as long as the cylinder radius is greater than zero. 

The proof is based on a novel semi-discrete, operator splitting 
numerical scheme, known as the kinematically coupled scheme, in- 
troduced in [32j to numerically solve the underlying FSI problems. 
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The backbone of the kinematically coupled scheme is the well-known 
Marchuk-Yanenko scheme, also known as the Lie splitting scheme. We 
effectively prove convergence of that numerical scheme to a solution 
of the corresponding FSI problem. 



1 Introduction 



We study the existence of a weak solution to a nonlinear moving bound- 
ary, unsteady, fluid-structure interaction (FSI) problem between an incom- 
pressible, viscous, Newtonian fluid, flowing through a cylindrical 2D domain, 
whose lateral boundary is modeled as a cylindrical Koiter shell. See Fig- 
ure [TJ Two Koiter shell models are considered: the linearly viscoelastic and 
the linearly elastic model. 

The fluid flow is driven by the time-dependent inlet and outlet dynamic 
pressure data. The fluid and structure are fully coupled via the kinematic 
and dynamic lateral boundary conditions describing continuity of velocity 
(the no-slip condition), and balance of contact forces at the fluid-structure 
interface. 



Thin Koiter shell 




Viscous, incompressible, Newtonian fluid 




Dynamic pressure inlet and outlet data 




Figure 1: Domain sketch. 



This work was motivated by a study of blood flow through medium-to- 
large human arteries. In medium-to-large human arteries, such as the aorta 
or coronary arteries, blood can be modeled as an incompressible, viscous, 
Newtonian fluid. Due to the high complexity and nonlinearity of the un- 
derlying FSI problem, various simplifications of the arterial wall mechanical 
properties have to be taken into account. The viscoelastic Koiter shell model, 
considered in this manuscript, and its many versions that appear in litera- 
ture, have been shown to be a good approximation of the mechanical prop- 
erties of arterial walls [9j [TO]. In fact, a version of the problem considered 
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in this manuscript has been used as a benchmark problem for FSI solvers 
in hemodynamics, see e.g., [27] . Therefore, studying its well-posedness has 
implications for both analysis and numerics. 

The development of existence theory for moving boundary, fluid-structure 
interaction problems, has become particularly active since the late 1990's. 
The first existence results were obtained for the cases in which the structure 
is completely immersed in the fluid, and the structure was considered to be 
either a rigid body, or described by a finite number of modal functions. See 
e.g., H7J US [2H |22l [221 123 IS], and the references therein. More recently, 
the coupling between the 2D or 3D Navier-Stokes equations and 2D or 3D 
linear elasticity on fixed domains, was considered for linear models in |24J, 
and for nonlinear models in [21 [31 36J. 

Concerning compliant (elastic or viscoelastic) structures, the first FSI ex- 
istence result, locally in time, was obtained in [I], where a strong solution for 
an interaction between an incompressible, viscous fluid in 2D and a ID vis- 
coelastic string was obtained, assuming periodic boundary conditions. This 
result was extended by Lequeurre in |37J, where the existence of a unique, 
local in time, strong solution for any data, and the existence of a global 
strong solution for small data, was proved in the case when the structure 
was modeled as a clamped viscoelastic beam. D. Coutand and S. Shkoller 
proved existence, locally in time, of a unique, regular solution for an inter- 
action between a viscous, incompressible fluid in 3D and a 3D structure, 
immersed in the fluid, where the structure was modeled by the equations 
of linear [18], or quasi- linear [19] elasticity. In the case when the structure 
(solid) is modeled by a linear wave equation, I. Kukavica and A. Tufahha 
proved the existence, locally in time, of a strong solution, assuming lower 
regularity for the initial data [35J. A fuid-structure interaction between a 
viscous, incompressible fluid in 3D, and 2D elastic shells was considered in 
[TBI IT2] where existence, locally in time, of the unique regular solution was 
proved. All the above mentioned existence results for strong solutions are lo- 
cal in time. We also mention that the works of Shkoller et al., and Kukavica 
at al. were obtained in the context of Lagrangian coordinates, which were 
used for both the structure and fluid problems. 

In the context of weak solutions, the following results have been obtained. 
Continuous dependence of weak solutions on initial data for a fluid structure 
interaction problem with a free boundary type coupling condition was stud- 
ied in [33]. Existence of a weak solution for a FSI problem between a 3D 
incompressible, viscous fluid and a 2D viscoelastic plate was considered by 
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Chambolle et al. in [TTJ, while Grandmont improved this result in [3T] to 
hold for a 2D elastic plate. In these works existence of a weak solution 
was proved for as long as the elastic boundary does not touch "the bottom" 
(rigid) portion of the fluid domain boundary. 

In the present manuscript we prove the existence of a weak solution to 
a FSI problem modeling the flow of an incompressible, viscous, Newtonian 
fluid flowing through a cylinder whose lateral wall is described by the lin- 
early viscoelastic, or by the linearly elastic Koiter shell equations. The fluid 
domain is two-dimensional, while the structure equations are in ID. The two 
existence results (the viscoelastic case and the elastic case) hold for as long 
as the compliant tube walls do not touch each other. The main novelty 
of this work is in the methodology of proof. The proof is based on a semi- 
discrete, operator splitting Lie scheme, which was used in [32] for a design of 
a stable, loosely coupled numerical scheme, called the kinematically coupled 
scheme (see also [7j). Therefore, in this work, we effectively prove that the 
kinematically coupled scheme converges to a weak solution of the underlying 
FSI problem. To the best of our knowledge, this is the first result of this kind 
in the area of nonlinear FSI problems. Semi-discretization is a well known 
method for proving existence of weak solutions to the Navier-Stokes equa- 
tions, see e.g. [17], Ch. III. 4. The Lie operator splitting scheme, also known 
as the Marchuk-Yanenko scheme, has been widely used in numerical compu- 
tations, see [30J and the references therein. Temam used a combination of 
these approaches in [16] to prove the existence of a solution of the nonlinear 
Carleman equation. The present manuscript represents the first use of this 
methodology in the area of nonlinear FSI problems. This method is robust 
in the sense that it can be applied to the viscoelastic case and to the elas- 
tic case, independently. Namely, our existence result in the case when the 
structure is purely elastic in not obtained in the limit, as the regularization 
provided by structural viscosity tends to zero. 

Another novelty of this work is in the the boundary conditions, which 
are not periodic, but are motivated by the blood flow application and are 
given by the prescribed inlet and outlet dynamic pressure data. This means 
that the Lagrangian framework for the treatment of the entire coupled FSI 
problem cannot be used in this context, and so we employed the Arbitrary 
Langrangian-Eulerian (ALE) method to deal with the motion of the fluid 
domain. 

Our proof is constructive, and its main steps have already been imple- 
mented in the design of several stable computational FSI schemes for the sim- 
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ulation of blood flow in human arteries [7J E2J [33] . The main steps in the 
proof include the ALE weak formulation and the time-discretization 
via Lie operator splitting. Solutions to the time-discretized problems de- 
fine a sequence of approximate solutions to the continuous time-dependent 
problem. The time-dependent FSI problem is discretized in time in such 
a way that at each time step, this multi-physics problem is split into two 
sub-problems: the fluid and the structure sub-problem. However, to achieve 
stability and convergence of the corresponding numerical scheme [32], the 
splitting had to be performed in a special way in which the fluid sub-problem 
includes structure inertia (and structural viscosity in the viscoelastic case) 
via a "Robin- type" boundary condition. See [5] for more details. The fact 
that structure inertia (and structural viscosity) are included implicitly in the 
fluid sub-problem, enabled us, in the present work, to get appropriate energy 
estimates for the approximate solutions, independently of the size of time 
discretization. Passing to the limit, as the size of the time step converges 
to zero, is achieved by the use of compactness arguments and a careful con- 
struction of the appropriate test functions associated with moving domains. 
The main difference between the viscoelastic and the elastic case is in the 
compactness argument. 

The main body of the manuscript is dedicated to the proof in the case 
when the structure is modeled as a linearly viscoelastic Koiter shell. In 
Section [8] we summarize the main steps of the proof in the case when the 
structure is modeled as a linearly elastic Koiter shell. 

2 Problem description 

We consider the flow of an incompressible, viscous fluid in a two-dimensional, 
symmetric cylinder (or channel) of reference length L, and reference width 
2R, see Figure [2j 

Without loss of generality, we consider only the upper half of the fluid 
domain supplemented by a symmetry boundary condition at the bottom 
boundary. Thus, the reference domain is (0, L) X (0, R) with the lateral (top) 
boundary given by (0, L) x {R}. 

We will be assuming that the lateral boundary of the cylinder is de- 
formable and that its location is not known a priori, but is fully coupled to 
the motion of the viscous, incompressible fluid occupying the fluid domain. 
Furthermore, it will be assumed that the lateral boundary is a thin, isotropic, 
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Figure 2: Domain sketch and notation. 



homogeneous structure, whose dynamics is modeled by the cylindrical, lin- 
early viscoelastic, or by the cylindrical linearly elastic Koiter shell equations. 
Additionally, for simplicity, we will be assuming that only the displacement 
in vertical (radial) direction is non-negligible. The (vertical) displacement 
from the reference configuration will be denoted by r)(t,z). See Figure [2] 
Models of this kind are common in blood flow applications [U H2J HI] , where 
the lateral boundary of the cylinder corresponds to arterial walls. 

The fluid domain, which depends on time and is not known a priori, will 
be denoted by 

n v (t) = {(z, r) G M 2 : z G (0, L), r G (0, R + T](t, x)}, 

and the corresponding lateral (top) boundary by 

T v (t) = {(z,r) G M 2 : r = R + rj{t, z), z G (0, L)}. 

The "bottom" (symmetry) boundary of the fluid domain will be denoted by 
Tb = (0, L) x {0}, while the inlet and outlet sections of the fluid domain 
boundary by r in = {0} x (0,R), T out = {L} x (0,R). See Figure § 

The fluid problem: We are interested in studying a dynamic pressure- 
driven flow through Q v (t) of an incompressible, viscous fluid modeled by the 
Naiver-Stokes equations: 

p^u + u-Vu) = V.,, j illfWiMe(0iT)i (1) 

where pf denotes the fluid density, u fluid velocity, p fluid pressure, a = 
—pi + 2/iD(u) is the Cauchy stress tensor of the fluid, p is the kinematic 
viscosity coefficient, and D(u) = |(Vu + V r u) is the symmetrized gradient 
of u. 
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At the inlet and outlet boundaries we prescribe zero tangential velocity 
and dynamic pressure p+ |w| 2 (see e.g. [16]): 



where Pin/out G Lf oc (0, oo) are given. Therefore the fluid flow is driven by a 
prescribed dynamic pressure drop, and the flow enters and leaves the fluid 
domain orthogonally to the inlet and outlet boundary. 

At the bottom boundary we prescribe the symmetry boundary condition: 



The structure problem, namely, the dynamics of the lateral bound- 
ary, is defined by the linearly viscoelastic cylindrical Koiter shell equations 
capturing radial displacement rj (for the purely elastic problem see Section [8]): 



Ps hd 2 t r] + Cov - C^rt + C 2 d 4 z v + Dodtr, - D&tijfi + D 2 d t d 4 z v = /• (4) 



Here, p s is the structure density, h is the structure thickness, and / is the 
force density in the radial (vertical) e r direction acting on the structure. 
The constants Cj and Dj > are the material constants describing structural 
elasticity and viscosity, respectively, which are given in terms of four material 
parameters: the Young's modulus of elasticity E, and the Poisson ratio er, 
and their viscoelastic couter-parts E v and a v (for a derivation of this model 
and the exact form of the coefficients, please see Appendix, and E]). The 
purely elastic case, i.e., Di = 0,i = 0, 1,2, will be considered in Section |8j 

We consider the dynamics of a clamped Koiter shell with the boundary 
conditions 



The coupling between the fluid and structure is defined by two sets of 
boundary conditions satisfied at the lateral boundary T v (t). They are the 
kinematic and dynamic lateral boundary conditions describing continuity of 
velocity (the no-slip condition), and continuity of normal stress, respectively. 
Written in Lagrangian framework, with z G (0, L), and for t G (0,T), they 
read: 

• The kinematic condition: 





u r = d r u z = 0, on T b . 



(3) 



?7(0) = 0*77(0) = r){L) = d zV (L) = 0. 



(d t r)(t,z),0) 



u(t,z,R + T)(t,z)), 



(5) 
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• The dynamic condition: 

f(t,z) = -J(t, z)(crn)\ {tiZ:R+r!(ttZ)) -e r . (6) 

Here, / = f(t, z) corresponds to the right hand-side of equation Q, and 
J(t, z) = \j\ + (d z r)(t, z)) 2 denotes the Jacobian of the transformation 
from Eulerian to Lagrangian coordinates. 

System ([T])-([6]) is supplemented with the following initial conditions: 

u(0,.)=uo, v(P,-)=m, d t 7 } (p,.)=v . (7) 

Additionally, we will be assuming that the initial data satisfies the fol- 
lowing compatibility conditions: 

u (z, R + r) (z)) = v (z)e r , z e (0, L), 

Vo (0) = 770 (L) = v (0) = v (L) = 0, (8) 
R + 7] (z)>0, ze[0,L]. 

Notice that the last condition requires that the initial displacement is such 
that the lateral boundary does not touch the bottom of the domain. This is 
an important condition which will be used at several places throughout this 
manuscript. 

In summary, we study the following fluid-structure interaction problem: 

Problem 2.1. Find u = (u z (t, z,r),u r (t, z,r)),p(t, z,r), and rj(t,z) such 
that 



p f (d t u + (u- V)u) = V-a 



V-u = 



mtt v (t), t e (0,T), 



(9) 



u 

p s hdfr] + C i] - Cidln + C 2 d z r] 
+D d t r] - Dxdtdln + D 2 d t d 4 zV 



d t ne r 



on (0,T) x (0,L), 



-Jan ■ e r , 



(10) 



u r = 0, 
d r u z = 0, 



on (0,T) x T b 



(11) 
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P + " lU j 2 _ ^/° Mt(t) ' ) on (0, T) x T in/out , (12) 
u(0, .) = u , ] 

77(0, .) = Vo, >aU = 0. (13) 
d t T](0, .) = v Q . ) 

This is a nonlinear, moving-boundary problem, which captures the full, 
two-way fluid-structure interaction coupling. The nonlinearity in the prob- 
lem is represented by the quadratic term in the fluid equations, and by the 
nonlinear coupling between the fluid and structure defined at the lateral 
boundary T v (t), which is one of the unknowns in the problem. 



2.1 The energy of the problem 



Problem (2.1) satisfies the following energy inequality: 

jE(t) + D{t) < C(P m (t), P out (t)), (14) 

where E(t) denotes the sum of the kinetic energy of the fluid and of the 
structure, and the elastic energy of the Koiter shell: 



P/ iiuii 2 + Psh \\a n \\i 

yll U llL2(Q„(t)) + -Y\\°tv\\mo,L) 

-L/^mmo M ^ no ^,11 11 o \ 



E(t) = ^l|u||i Wt) ) + ^ll^ 

+ ^(^II^^W) + C lW d ^\\h(0,L) + C 2\\ d zV\\h(0,L))i 

the term D{t) captures dissipation due to structural and fluid viscosity 



D(t) =^||D(u) „(*«)H-^o||^^IIf^ C o.i>^^iIl^^^lli^Co,x,5H-^2ll^t^^ll^co,i>» 

(16) 

and C (Pi n (t) , P out (t))) is a constant which depends only on the inlet and 
outlet pressure data, which are both functions of time. 

To show that (14) holds, we first multiply equation by u, integrate 



over ft„(t) } and formally integrate by parts to obtain: 
/ p / (9 i u-u+(u-V)u-u)+2 / u / |Du| 2 - / (-pI+2//D(u))n(i)-u = 0. 

To deal with the inertia term we first recall that Q^t) is moving in time and 
that the velocity of the lateral boundary is given by u|m)- The transport 
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theorem applied to the first term on the left hand-side of the above equation 
then gives: 

d t u-u=-— [ |u| 2 / |u| 2 u-n(t). 



The second term on the left hand side can can be rewritten by using inte- 
gration by parts, and the divergence-free condition, to obtain: 



/ (u-V)u u = J/ |u| 2 u-n(t) = ^( f |u| 2 u-n(t) 

/ |u| 2 m z .) 

•J r 'out 



\u\ 2 u. 



These two terms added together give 

/ 3jU-u+ / (u-V)u-u = ~ / |u| 2 -^ / |u| 2 u 2 + ^ / \u\ 2 u z . 
Jn v (t) Jn v (t) Z( ttJn v (t) z Jr in z Jr out 

(17) 

To deal with the boundary integral over dQ n (t), we first notice that on 
Tin/out the boundary condition ^ implies u r = 0. Combined with the 
divergence-free condition we obtain d z u z = —d r u r = 0. Now, using the 
fact that the normal to T in / out is n = (=)=1, 0) we get: 

I (-pI + 2//D(u))n-u= / P m u z - I P out u z . (18) 

in/ out JT'in •> T 'out 

In a similar way, using the symmetry boundary conditions ([3]), we get: 

[-pi + 2//D(u))n ■ u = 0. 



What is left is to calculate the remaining boundary integral over r v (t). 
For this we consider the Koiter shell equation Q, multiply it by d t r], and 
integrate by parts to obtain 

f L fa P shd \\A 112 

+\f t (0)IMl! W ) + Ci||^||| 2(0)i) + C 2 ||^||| 2(0)i) ) (19) 



+D \\d t T]\\l 2i0jL) + ^111^^1^2(0^) + ^211^771^2(0^). 
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By enforcing the dynamic coupling condition ([6]) we obtain 



<rn(t) ■ u — — Jan ■ u 
r„(t) Jo 



fdtV- 



(20) 



Finally, by combining (20) with (19), and by adding the remaining con 



tributions to the energy of the FSI problem calculated in equations (17) and 

1 d 



(18), one obtains the following energy equality: 



P±d_ 
2 dt 



I |2 , P* h d \\a I 

|u| + ^Tdt lldtr]l 



l- n i.l i rl^i I |Du| z + ~(C ||77||| 2(0ii) 



fi„(t) 



+Ci\\dzr)\\l*(o,L) + C 2\\d 2 z v\\l^o,L)) + A)||9t*7||£a ( o, 



Di\\d t d z T]\\l2 {QjL) + D 2 \\d t dlri\\ 2 L 2 



(0,L) 



in j out 



(t) 



(21) 



in/ out 



By using the trace inequality and Korn inequality one can estimate: 



in j out 



c 

U z \ < C\P in / out \ ||u||^i(^(t)) < 



n,/ out 



out 



2 +^-\\V(u)\\ 2 L2(av(t) . 



By choosing e such that f</iwe get the energy inequality (|14|). 



3 The ALE formulation and Lie splitting 

3.1 First order ALE formulation 



To prove the existence of a weak solution to Problem 2.1 it is convenient to 



map Problem 2.1 onto a fixed domain f2. In our approach we choose f2 to 
be the reference domain Q = (0, L) x (0, 1). We follow the approach typical 
of numerical methods for fluid-structure interaction problems and map our 
fluid domain Q(t) onto Q by using an Arbitrary Lagrangian-Eulerian (ALE) 
mapping [7J [32j [231 1131 HI] • We remark here that in our problem it is not 
convenient to use the Lagrangian formulation of the fluid sub-problem, as is 
done in e.g., P32 Q21 ES], since, in our problem, the fluid domain consists of 
a fixed, control volume of a cylinder, which does not follow Largangian flow. 
We begin by defining a family of ALE mappings A n parameterized by 77: 



A v (t) : Q ->• fi„(t), A v {t)(z,r) :-- 



(R + r}{t,z))r 



[~z,f)en, (22) 
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where (z, r) denote the coordinates in the reference domain Q = (0, L) x (0, 1). 
Mapping A n (t) is a bijection, and its Jacobian is given by 

\deWA ri (t)\ = \R + r](t,z)\. (23) 

Composite functions with the ALE mapping will be denoted by 

uP{t,.)=u(t,.)oA n (t) and p^t,.) =p(t,.)oA v (t). 

The derivatives of composite functions satisfy: 

d t u = dtvP - (w" ■ V"K, Vu = V'u", 

where the ALE domain velocity, w v , and the transformed gradient, V, are 
given by: 

' Oy — r— df » 



w r ' = d t rjfe r , V 



1 



/ 



Note that 



V"v = Vv(V^)- 1 . 
The following notation will also be useful: 

= -p*I + 2/iDV), D^u") = -(Vu" + (V) T u r ' 



(24) 



(25) 



We are now ready to rewrite Problem |2.1| in the ALE formulation. However, 
before we do that, we will make one more important step in our strategy to 



prove the existence of a weak solution to Problem 2.1 Namely, as mentioned 
in the Introduction, we would like to "solve" the coupled FSI problem by 
approximating the problem using time-discretization via operator splitting, 
and then prove that the solution to the semi-discrete problem converges to 
a weak solution of the continuous problem, as the time-discretization step 
tends to zero. To perform time discretization via operator splitting, which 
will be described in the next section, we need to write our FSI problem 
as a first-order system in time. This will be done by replacing the second- 
order time-derivative of rj, with the first-order time-derivative of the structure 
velocity. To do this, we further notice that in the coupled FSI problem, the 
kinematic coupling condition (J5f implies that the structure velocity is equal 
to the normal trace of the fluid velocity on T v (t). Thus, we will introduce a 
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new variable, v, to denote this trace, and will replace d t r) by v everywhere 
in the structure equation. This has deep consequences both for the existence 
proof presented in this manuscript, as well as for the proof of stability of the 
underlying numerical scheme, presented in [8], as it enforces the kinematic 
coupling condition implicitly in all the steps of the scheme. 

Thus, Problem 2.1 can be reformulated in the ALE framework, on the ref- 
erence domain Q, and written as a first-order system in time, in the following 
way: 

Problem 3.1. Find u(t, z,r),p(t, z,r),r](t, z) , andv(t,z) such that 



p f (d t u+ ((u- w^) • V)u) = V-o* 



drltz 



P + P 4\u\ 2 

U r 



V -u 

= 



in (0,T) x tt, 







' on(0,T)xr 6 



Pin/out (t) j 



on (0,T) x T in/out , 



(26) 
(27) 
(28) 



u 

d t r] 

p s hd t v + C rj - Cid 2 z r] + C 2 d*ri 
+D v - D x d 2 z v + D&v 

u(0, .) =u o ,r}(0,.) =r]o,v(0, .) =v , at t = 0. (30) 

Here, we have dropped the superscript rj in u v for easier reading. 

We are now ready to define the time discretization by operator splitting. 
The underlying multi-physics problem will be split into the fluid and struc- 
ture sub-problems, following the different "physics" in the problem, but the 
splitting will be performed in a particularly careful manner, so that the re- 
sulting problem defines a scheme which converges to a weak solution of the 
continuous problem (and provides a numerical scheme which is uncondition- 
ally stable). 



on (0,T) x (0,L), (29) 



- Jan ■ e r 
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3.2 The operator splitting scheme 

We use the Lie splitting, also known as the Marchuk-Yanenko splitting 
scheme. The splitting can be summarized as follows. Let N eN, At = T/N 
and t n = nAt. Consider the following initial-value problem: 

^ + A0 = O in(0,T), 0(0) = O , 

where A is an operator defined on a Hilbert space, and A can be written as 
A = Ai + A 2 . Set 0° = 0o, an d, for n = 0, . . . , iV — 1 and i = 1,2, compute 
n +2 by solving 

at . . > m {t n ,t n+ i), 



4>i(t r 



and then set n+ 2 = 0j(t ra+1 ), for i = 1,2. It can be shown that this method 
is first-order accurate in time, see e.g., [30] . 



We apply this approach to split Problem |3.1| in the fluid and structure 
sub-problems. During this procedure the structure equation (29) will be 
split into the viscous part, i.e., the part involving the normal trace of the 
fluid velocity on T v (t), v, and the purely elastic part. The viscous part of 
the structure problem will be used as a boundary condition for the fluid 
sub-problem, while the elastic part of the structure problem will be solved 



separately. More precisely, we define the splitting of Problem 3.1 in the 
following way: 

Problem Al: The structure elastodynamics problem. In this 
step we solve the elastodynamics problem for the location of the deformable 
boundary by involving only the elastic energy of the structure. The motion 
of the structure is driven by the initial velocity, which is equal to the trace of 
the fluid velocity on the lateral boundary, taken from the previous step. The 
fluid velocity u remains unchanged in this step. More precisely, the problem 
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reads: Given (u n , r] n , v n ) from the previous time step, find (u, v, rj) such that: 
d t u = 0, in (t n ,t n+1 ) x CI, 

p s hd t v + C r] - C-ffi) + C 2 djr] = on (t n , t n+1 ) x (0, L), 
< d t 7] = v on (t n ,t n+1 ) x (0,L), (31) 
77(0) = 0^(0) = 77(L) = g> 2?7 (L) = 0, 

k u(t n ) = u n , n(t n ) = V n , v(t n ) = v n . 

Then set u™+§ = u(t n+l ), rf*h = n(t n+1 ), v n ^ = v(t n+1 ). 

Problem A2: The fluid problem. In this step we solve the Navier- 
Stokes equations coupled with structure inertia and viscoelastic energy of the 
structure, through a "Robin-type" boundary condition on T. The kinematic 
coupling condition is implicitly satisfied. The structure displacement remains 
unchanged. With a slight abuse of notation, the problem can be written as 
follows: Find (u,v,rj) such that: 



d t r) = on (t n , t n+1 ) x (0, L), 

„(*» + «».- w^).vn«) = y»"^- 1 in(t „, t „ +l)x(i , 

p,hS,v + D a v - Did'(v + D 2 div = - Jan ■ e r 1 on (tnj(n+l) x 

d r u z = o ) ™(*n,Wi)xr 6 , 

P+^\u\ 2 = Pin/outit) 1 , ,WF. , 

M _ g ? Oli ^t n , tri+lj x 1 in/out: 

withu(t„,.) = u" + i, 77(t n ,.)=r7" + i, t;(t n ,.)=^+i (32) 

Tien set u" +1 = u(t n+1 ), 77™ +1 = r)(t n+1 ), v n+l = v(t n+1 ). 

Notice that, since in this step 77 does not change, this problem is linear. 
Furthermore, it can be viewed as a stationary Navier-Stokes-like problem on 
a fixed domain, coupled with the viscoelastic part of the structure equation 
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through a Robin-type boundary condition. In numerical simulations, one can 
use the ALE transformation A v n to "transform" the problem back to domain 
Q v n and solve it there, thereby avoiding the un-necessary calculation of the 
transformed gradient V™. The ALE velocity is the only extra term that 
needs to be included with that approach. See, e.g., [7j for more details. For 
the purposes of our proof, we will, however, remain in the fixed, reference 
domain Q. 

It is important to notice that in Problem A2, the problem is "linearized" 
around the previous location of the boundary, i.e., we work with the domain 
determined by r/ n , and not by ^ n+1 / 2 . This is in direct relation with the 
implementation of the numerical scheme studied in [7J [SJ. However, we also 
notice that ALE velocity, w n+ 2, is taken from the just calculated Problem 
Al! This choice is crucial for obtaining a semi-discrete version of an energy 
inequality, discussed in Section [5] 

In the remainder of this paper we use the splitting scheme described above 
to define approximate solutions of Problem 3.1| (or equivalently Problem 2.1 ) 
and show that the approximate solutions converge to a weak solution, as 
At ->■ 0. 



4 Weak solutions 

4.1 Notation and function spaces 



To define weak solutions of the moving-bounday problem |3.1 we first intro- 



duce some notation which will simplify the subsequent analysis. We begin 
by introducing the following bilinear forms associated with the elastic and 
viscoelastic energy of the Koiter shell: 

a s (v, = r (C rpl> + C&rjdrf + C 2 d 2 zV d 2 j) , (33) 
Jo 

a'sfa tfO = I { D oV<P + D 1 d z r j d z iP + D 2 d 2 V d 2 ^) . (34) 
Jo 

Furthermore, we will be using b to denote the following trilinear form cor- 
responding to the (symmetrized) nonlinear term in the Navier-Stokes equa- 
tions: 

b(t,u, v,w) = - / (u • V)v ■ w / (u • V)w • v. (35) 

2 Jn„(t) 2 Jn ft) 
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Finally, we define a linear functional which associates the inlet and outlet 
dynamic pressure boundary data to a test function v in the following way: 

( F (t), v)r in/out = Pin(t) f v z - P out (t) I v z . 



To define a weak solution to Problem |3.1| we introduce the necessary 
function spaces. For the fluid velocity we will need the classical function 
space 

V F {t)= {u=(u z ,u r )eH 1 (n v (t)) 2 :V-u = 0, 

u z = on T(t), u r = on Q v (t) \ T(t)}. 1 ' 

The function space associated with weak solutions of the Koiter shell is given 
by 

V s = H 2 (0,L). (37) 

Motivated by the energy inequality we also define the corresponding evolution 
spaces for the fluid and structure sub-problems, respectively: 

W F (0, T) = L°°(0, T; L 2 (n v (t)) D L 2 (0, T; V F (t)) (38) 

W s (0,T) = W 1 ' oo (0,T;L 2 (0,L))n J ff 1 (0,T;Vs). (39) 

The solution space for the coupled fluid-structure interaction problem must 
involve the kinematic coupling condition. Thus, we define 

W(0,T) = {(11,77) e W F (0,T) x W s (0,T) : u(t, i2 + »j(f, z)) = d tV (t, z)e r }. 

(40) 

The corresponding test space will be denoted by 

Q(0, T) = {(q, ^) G ^([0, T); V F x V s ) : q(t, 2, i? + *)) = z)e r }. 

(41) 

4.2 Weak solution on the moving domain 

We are now in a position to define weak solutions of our moving-boundary 
problem, defined on the moving domain £l v (t). 
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Definition 4.1. We say that (11,77) G W(0, T) is a weak solution of Prob- 
lem 2.1 if for every (q,ip) € <2(0,T) the following equality holds: 

r-T 



Pf{- f [ u d t q+ [ 6(f,u,u,q))+2/z / / D(u) : D(q) 

J0 JSl,,(i) JO JO J!l,,(i) 

-^(1 (dtvYiJ - Psh [ [ d t r]d^+ I [asi^ + a'sidtr]^)) 



Jo Jo Jo 



= [ (F(t),q)r in/out + Pf [ u -q(0)+p s h[ ^(0). 
Jo Jn, 10 Jo 

(42) 

In deriving the weak formulation we used integration by parts in a classical 
way, and the following equalities which hold for smooth functions: 

/ (u-V)uq = \! (u-V)uq-J/ (u • V)q • u 

Jn v (t) A Jn v (t) 1 Jn v (t) 



out j in 

T pL 



/ / d t u ■ q = - / / u d t q- u • q(0) - / / 

Jo Jn v (t) Jo Jn v (t) Jn m) Jo Jo 



(d t rj) ip- 



4.3 Weak solution on a fixed, reference domain 

Since most of our analysis will be performed on the problem defined on the 
fixed, reference domain Q, we rewrite the above definition in terms of Q using 
the ALE mapping A v (t) defined in (22). For this purpose, we introduce the 
notation for the transformed trilinear functional b v , and the function spaces 
for the composit, transformed functions defined on the fixed domain Q. 
The transformed trilinear form W is defined as: 

6 V, u, q) = i / (R + V ) (((u - w') ■ V> • q - ((u - w') ■ V)q • u) , (43) 

where R + rj is the Jacobian of the ALE mapping, calculated in (23). Notice 
that we have included the ALE domain velocity w v into W. 
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It is important to point out that the transformed fluid velocity u v is not 
divergence-free anymore. Rather, it satisfies the transformed divergence-free 
condition V • \i n = 0. Therefore we need to redefine the function spaces for 
the fluid velocity by introducing 

V F = {u = (u z ,u r ) E H\n) 2 : V 7 • u = 0, u z = on T, u r = on \ T}. 

The function spaces W F (0,T) and W v (0,T) are defined the same as before, 
but with Vp instead VV(t). More precisely: 

W£(0,T) = L°°(0,T;L 2 (fi)nL 2 (0,T;V^(t)), (44) 

W(0,T) = {(11,77) e W^(0,T) x Ws(0,T) : u(t,z,l) = <9 t r/(t, z)e r }. (45) 
The corresponding test space is defined by 

Q?{0,T) = {(q,V) e ^([0,1)^; x V s ) : q{t,z,l) = i>(t,z)e r }. (46) 

Definition 4.2. H^e say (u, rj) e VV(0, T) zs a u>ea£; solution of Prob- 



lem 3.1 defined on the reference domain VL, if for every (c[,ip) G Q v (0,T) the 
ty holds: 

Pf(~ f ! (R + v)u-d t q+ /V( U ,u,q)) 
Jo Jn Jo 

2/i f [ (R + V )W(u) :D"(q)-^ f [ (d tV )u ■ q 
- Ps h / d tV d t iP+ (a s (v^) + a' s (d tV ,i/j)) (47) 



following equality holds. 



jo 



R / [Pin{t) / (g 2 )| 2 =0 ~~ Pout 
Jo Jo Jo 

u • q(0) + p s h / voip(0). 
Jo 



+Pf 

To see that this is consistent with the weak solution defined in Defini- 



tion 4.1, we present the main steps in the transformation of the first integral 
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on the left hand-side in (42 ), responsible for the fluid kinetic energy. Namely, 
we formally calculate: 

- / u • d t q = - I (R + 77)11" ■ (d t q - (w" ■ V)q) = - I (R + 77)11" ■ <9 t q 
JUr, Jn Jn 



1 



+ / ( R + V)(™ V ■ V)q • u" + - / (R + 77) (w" • V)q • u". 
^ ^ Jn 

In the last integral on the right hand-side we use the definition of w" and of 

V, given in (24), to obtain 

( J R + r])(w' ? ■ V)q- u" = / ^f3 f q-u". 
! Jn 

Using integration by parts with respect to r, keeping in mind that 77 does not 
depend on r, we obtain 



u • d t q 



(r + r/K ■ ($q - (W • V) q ) = - (R+ r/K ■ 9tq 



+ / OR+r/Xw^V^q-u^- 



( J R+/7)(w ?? -V ,7 )u ? '-q- 



^u^-q-l-; 



By using this identity in (42), and by recalling the definitions for b and W, 
we obtain exactly the weak form (47). 

In the remainder of this manuscript we will be working on the fluid- 
structure interaction problem defined on the fixed domain Q, satisfying the 
weak formulation presented in Definition 4.2 For brevity of notation, since 
no confusion is possible, we omit the superscript "tilde" which is used to 
denote the coordinates of points in Q. 



5 Approximate solutions 

In this section we use the Lie operator splitting scheme and semi-discretization 



to define a sequence of approximations of a weak solution to Problem 3.1 



Each of the sub-problems defined by the Lie splitting in Section 3.2| (Problem 



Al and Problem A2), will be discretized in time using the Backward Euler 
scheme. This approach defines a time step, which will be denoted by At, and 
a number of time sub-intervals iV e N, so that 



(0,T) = ut- 1 (r,r +1 ), t n = nAt, n = 0,...,N-l. 
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For every subdivision containing JVeN sub-intervals, we recursively define 
the vector of unknown approximate solutions 



N 



( u «+i \ 



V N ' 
n+l 

\ Vn ' 



,71 



0,1, 



N-h 



1,2, 



(48) 



/ 



where i — 1,2 denotes the solution of sub-problem Al or A2, respectively. 
The initial condition will be denoted by 



X° 



The semi-discretization and the splitting of the problem will be performed 




in such a way that the discrete version of the energy inequality (14) is pre 



served at every time step. This is a crucial ingredient for the existence proof. 
We define the semi-discrete versions of the kinetic and elastic energy, 



originally defined in (15), and of dissipation, originally defined in (16), by 
the following: 



E n+±2 

N 



n-l+i\i n +f |2 



U 



A? 



+ p s h\\v 



2 II 2 
N llL 2 (0,L) 



"Co||% 2 1| |a (0,2,) + Ci\\d z 7) N 2 ||i2( 0jI/ ) + C 2 \\d 2 z rj N \\ L - l(lL] j- 



(49) 



jjn+l 



At(ji J {R + v 

4-D llr) 2 7! n+1 ll 2 



D r ' n (u n N +1 )\ 2 



n iu, n+1 n 2 



n \\f) <N n+1 ii 2 

L>l\\O z V N || L 2 (0ii) 



z v N \\L 2 (0,L)h 



n 



0,...,N 



1, % 



0,1. 



(50) 

Throughout the rest of this section, we fix the time step At, i.e., we keep 
iV G N fixed, and study the semi-discretized sub-problems defined by the 
Lie splitting. To simplify notation, we will omit the subscript N and write 



(ll 2 , V 2 , Tj 2 



. r , n+5 n+£ n+f- 

instead of {u N ,v N ,7} N 



5.1 Semi-discretization of Problem Al 

We write a semi-discrete version of Problem Al (Structure Elastodynamics) , 



defined by the Lie splitting in (31). In this step u does not change, and so 



u 



u . 
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We define (v n+ 2 , r] n+ 2 ) e H$(0,L) x Hq(0,L) as a solution of the following 
problem, written in weak form: 



jf ^At r <f> = £ v n+1 H, 0GL 2 (O,L), 
— ^ + a s (r/ n+ 5,^) = 0, ^eH 2 (0,L). 



(51) 



The first equation is a weak form of the semi-discretized kinematic coupling 
condition, while the second equation corresponds to a weak form of the semi- 
discretized elastodynamics equation. 



Proposition 5.1. For each fixed At > 0, problem (51) has a unique solution 
(v n+ *,7} n+ *) G H 2 (0,L) x # 2 (0,L). 

Proof. The proof is a direct consequence of the Lax-Milgram Lemma applied 
to the weak form 

V n+ h + (At) 2 a s (v n+ K^)= [ L (Atv n + V n )if>, ^etf 2 (0,L), 



which is obtained after elimination of v n+ 2 in the second equation, by using 
the kinematic coupling condition given by the first equation. □ 



Proposition 5.2. For each fixed At > 0, solution of problem (51) satisfies 
the following discrete energy equality: 

El +l + \ -v n \\ 2 + C \\v n ^ ~ V n \\ 2 (52) 

+Ci||0,fa B+ i - rf )|| 2 + C 2 \\d 2 M +l ~ V n )\\ 2 ) = K, 

where the kinetic energy E 1 ^ is defined in (|49|. 



Proof. From the first equation in (51) we immediately get 

Therefore we can take t> n+ i as a test function in the second equation in (51 ). 
We replace the test function if) by v n+ 2 in the first term on the left hand-side, 
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and replace ip by {j] n+ ^ — r] n )/ At in the bilinear form as- We then use the 
algebraic identity (a — b) ■ a = |(|a| 2 + \a — b\ 2 — \b\ 2 ) to deal with the terms 
^n+1/2 _ v n^ v n+i/2 an j ^n+i/2 _ _ After multiplying the entire 



equation by At, the second equation in (51) can be written as 



p s ^(||?;" + *|| 2 + _ ^|| 2 ) + a^^+^^^+s) + 0^(77^+5 - 77^,^+5 - 77") 

= p s h\\v n \\ 2 + a s (v n ,V n )- 

We then recall that u n+ 5 = u ra in this sub-problem, and so we can add 
Pf J n (l + r] n )u n+1 / 2 on the left hand-side, and pj J Q (1 + r) n )u n on the right 
hand-side of the equation, to obtain exactly the energy equality (52). □ 



5.2 Semi-discretization of Problem A2 

We write a semi-discrete version of Problem A2 (The Fluid Problem), defined 



by the Lie splitting in (32). In this step 77 does not change, and so 



Define (u n+1 , v n+1 ) G Vjf x H$(0, L) by requiring that for each (q, ip) e V F x 
Hq(0, L) such that q|r = ^e r , the following weak formulation of problem (32) 
holds: 



Pf (R + rf) 
Jn 



u n+l _ u n+| 



1 



At 



• q 



(u n -v n+ 2re r ) ■ V< 



u n+1 • q 



q-u !l+1 ) + y J v n+1 *u n+1 ■ q 



[u n -v n+ ire r ) ■ V v 

+2/i/ n ( J R + 77 n )D" n (u):D" n (q) 



+p s h 



L v n+l _ v n+\ 



I /„,ra+l 



At 



tj; + a' s (v 



i>) = R(P? n [\q z )\ z= o ~ Pit I (q z )\ z =L): 
Jo Jo 



with V" • u n+1 = 0, u™ +1 = v n+1 e r 



where P? n/out = — / P in/out (t)dt. 

^ JnAt 



(53) 
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Proposition 5.3. Let At > ; and assume that if 1 are such that R + if 1 > 
i?min > 0, n = 0,...,N. Then, the fluid sub-problem defined by (53) has a 
unique weak solution (u n+1 ,v n+1 ) 6 V n F x Hq(0,L). 

Proof. The proof is again a consequence of the Lax-Milgram Lemma. More 
precisely, donote by U the Hilbert space 



U = {(u,v) E Vf x H 2 (0,L) : u, r = ve z }, 



(54) 



and define the bilinear form associated with problem (53): 

At 



a((u,u),(q,V)) := Pf / (R + rf) u ■ q 



At 

T 



(u n - v n+ ire r ) ■ V" 



u 



q ■ u 



u • q 



+ At 



Pf 



v n+ *u ■ q + At2fi / (i2 + )D" (u) : D" (q) 

+ p s /i / vip + Ata' s (v,ip), (u, f), (q, ■0) G U. 
Jo 

We need to prove that this bilinear form a is coercive and continuous on U. 
To see that a is coercive, we write 

f At i 

a({u,v),{u,v)) = p f / {R + v n + -7T v n+1 i)\u\ 2 + Ps h 



V)). 



+At(2p [ (R + V n )\W n (u)\ 2 + a' s (v, 
Jn 



Coercivity follows immediately after recalling that r] n are such that R + r] n > 
Rmm > 0, which implies that R+7] n + fv n+ ^ = R+j(r] n +7] n+1 2) > R min > 0. 
Before we prove continuity notice that from (24) we have: 

llV^ull^n) < C\\r] n \\ H 2 {0jL) \\u\\ H i {n) . 

Therefore, by applying the generalized Holder inequality and the continuous 
embedding of H 1 into L 4 , we obtain 

o((u,f), (q,V)) < C (p/||u|| £ 2(n)||q||L2(n) + p s h\\v\\ L 2 ( ^ L) \\iI)\\ l *{o,l) 



+At\\rj n \\ H 2 {QjL) {\\xi n \\ H i m + \\v n+1 2\\ H i (0jL) ) ||u|| ff i(n)||q||Hi(n) 
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+ At/i||r/ n ||^ 2(0iL) ||u||Hi(n)||q||//i(f7) + At|M|ff2 ( o,L)|M|ff2 ( o iL )J . 

This shows that a is continuous. The Lax-Milgram lemma now implies the 
existence of a unique solution (u n+1 , v n+1 ) of problem (53). □ 



Proposition 5.4. For each fixed At > 0, solution of problem (53) satisfies 
the following discrete energy inequality: 



E n N +l + Pl I ( J R + r? ™)|u" +1 -u n | 2 + ^||t; n+1 -^ + -" 2 



2 V" ' ■' " 1 ' 2 (55) 



where the kinetic energy E^ and dissipation are defined in (49) and (50), 
and the constant C depends only on the parameters in the problem, and not 
on At (or N). 



Proof. We begin by focusing on the weak formulation (53 ) in which we replace 
the test functions q by u n+1 and ip by v n+1 . We multiply the resulting 
equation by At, and notice that the first term on the right hand-side is given 
by 

|^ + rf)|u^| 2 . 

This is the term that contributes to the discrete kinetic energy at the time 
step n + 1, but it does not have the correct form, since the discrete kinetic 
energy at n + 1 is given in terms of the structure location at n + 1, and not 
at n, namely, the discrete kinetic energy at n + 1 involves 



^ I (R + V n+1 )\u n+1 ^ 



To get around this difficulty it is crucial that the advection term is present 
in the fluid sub-problem. The advection term is responsible for the presence 
of the integral 

P -L f Atv n+ hu n+1 \ 2 
2 Jn 11 

which can be re-written by noticing that Atv n+ ^ := (r/ ra+1 / 2 — rf 1 ) which is 
equal to (rf +l — rf) since, in this sub-problem r] n+1 = ^ n+1 / 2 . This implies 



(R + 7] n )\u n+1 \ 2 + Atv n+1 2\u n+1 \ 2 ^j = — f (R + rj n+1 )\u n+1 ^ 



25 



Thus, these two terms combined provide the discrete kinetic energy at the 
time step n+1. It is interesting to notice how the nonlinearity of the coupling 
at the deformed boundary requires the presence of nonlinear advection in 
order for the discrete kinetic energy of the fluid sub-problem to be decreasing 
in time, and to thus satisfy the desired energy estimate. 

To complete the proof one simply uses the algebraic identity (a — b) ■ a 
U\a\ 2 + \a-b\ 2 -\b\ 2 ) in the same way as in the proof of Proposition 



5.2 



□ 



We pause for a second, and summarize what we have accomplished so 
far. For a given At > we divided the time interval (0,T) into N = T/At 
sub-intervals (t n ,t n+1 ),n = 0,...,N — 1. On each sub-interval (t n ,t n+1 ) we 
"solved" the coupled FSI problem by applying the Lie splitting scheme. First 
we solved for the structure position (Problem Al) and then for the fluid flow 
(Problem A2). We have just shown that each sub-problem has a unique so- 
lution, provided that R + rf 1 > R m i n > 0, n = 0, N, and that its solution 
satisfies an energy estimate. When combined, the two energy estimates pro- 



vide a discrete version of the energy estimate (14). Thus, for each At we 



have a time-marching, splitting scheme which defines an approximate solu- 



tion on (0,T) of our main FSI problem defined in Problem 3.1, and is such 
that for each At the approximate FSI solution satisfies a discrete version of 
the energy estimate for the continuous problem. 

What we would like to ultimately show is that, as At — > 0, the sequence 
of solutions parameterized by iV (or At), converges to a weak solution of 



Problem |3.1| Furthermore, we also need to show that R + rf 1 > R m m > is 
satisfied for each n — 0, N — 1. In order to obtain this result, it is crucial 
to show that the discrete energy of the approximate FSI solutions defined for 
each At, is uniformly bounded, independently of At (or N). This result is 
obtained by the following Lemma. 

Lemma 5.1. (The uniform energy estimates) Let At > and N = 

T/At > 0. Furthermore, let E^ 2 ,E r ^ 1 , and D^ N be the kinetic energy and 



dissipation given by (49) and (50), respectively. 

There exists a constant C > independent of At ( and N ), which depends 
only on the parameters in the problem, on the kinetic energy of the initial 
data Eq, and on the energy norm of the inlet and outlet data \\Pi n /out\\ 2 i2^ T y 
such that the following estimates hold: 

1. E n ^ < C,E n N +1 < C, for alln = 0,...,N - 1, 
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*■ £f=i K < c, 



N-l 



3. V ( f(R+ rf)\vi n+1 - iff + \\v n+1 

n=0 ^ 



- l) n+ l\\ 2 

V \\L 2 (0,L) 



N-l 



l £ ((Coll^ 1 - r/"||i 2(0iL) + CxUW +1 - ^)|| 2 i2(0>L) 

+C 2 ||9 2 (^-^)||i 2(0iL) ) <C. 



n=0 



/fi/oci, C = Eq + C (\\Pin\\ i 2 (o,T) + ll-fo"tlli2(o,T)J > w here C is the constant 
from (55), which depends only on the parameters in the problem. 



Proof. We begin by adding the energy estimates (52) and (55) to obtain 



E n N +l + D n N +1 + - (p f [ (R + r) n )\u n+1 - u n \ 2 + p s h\\v n+1 - v n+ ^\ 
2 v Jn 



2 

L 2 (0,L)~ 



+ Ps h\\v n+ -> - v n \\l^ L) + Coh n+ -> ~ V n \\l H o,L) + CiWdiirr** - rf)||! W) + 

+C 2 ||a z 2 (^-^)||i 2(0ii) ) <E N +CAt((Pp n f+(P^ t ) 2 ), n = 0,...,N-l. 

Then we calculate the sum, on both sides, and cancel the same terms in the 
kinetic energy that appear on both sides of the inequality to obtain 



N-l N-l 



E N N + Y. D N +l + \Y.{pf S (R + V n )\" n+1 -U n \ 2 +Psh\\v 



n+l_ n+|||2 , 

w IIl 2 (o,l)+ 



n=0 n=0 

+MlK +f - v n \\l^ L) + Coll^ +i " ^IIW,L) + Cx||a.(^ - + 

JV-1 

+C 2 \\d 2 M^ - V n )\\h( ,L)) <E + CAtJ2((P- n ) 2 + (Pit?)- 



n=0 

To estimate the term involving the inlet and outlet pressure we recall that on 
every sub-interval (t n ,t n+l ) the pressure data is approximated by a constant 
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which is equal to the average value of the pressure over that time interval. 
Therefore, we have, after using Holder's inequality: 



N-l N-l 



n=0 




f(n+l)Ai 



nAt 



By using the pressure estimate to bound the right hand-side in the above 
energy estimate, we have obtained all the statements in the Lemma, with 
the constant C given by C = E + C\\P in / ou t\\ 2 L 2^ T y 

Notice that Statement 1 can be obtained in the same way by summing 
from to n — 1, for each n, instead of from to N — 1. □ 

We will use this Lemma in the next section to show convergence of ap- 
proximate solutions. 



6 Convergence of approximate solutions 



We define approximate solutions of Problem 3.1 on (0, T) to be the functions 
which are piece-wise constant on each sub-interval ((n — 1) At, nAt], n = 
1 ... AT of (0, T), such that for t E ((n - 1) At, nAt], n — l...N, 

n-l 

Ujv(t, .) = U^, T] N (t, .) = 7&, V N (t, .) = Vjf, v* N (t, .) = v N 2 . (56) 

See Figure 3. Notice that functions v* N = v 1 ^ 1//2 are determined by Step Al 
(the elastodynamics sub-problem), while functions vn = v 7 ^ are determined 
by Step A2 (the fluid sub-problem). As a consequence, functions vn are 
equal to the normal trace of the fluid velocity on T, i.e., u N = V]ye r . This 
is not necessarily the case for the functions v* N . However, we will show later 
that the difference between the two sequences converges to zero in L 2 . 



Using Lemma 5A we now show that these sequences are uniformly bounded 
in the appropriate solution spaces. 

We begin by showing that (r] N ) NeN is uniformly bounded in L°°(0, T; Hq(0, L)), 
and that there exists a T > for which R + rffr > holds independently of 
A^ and n. This implies, among other things, that our approximate solutions 
are, indeed, well-defined on a non-zero time interval (0, T). 

Proposition 6.1. Sequence (^^eN is uniformly bounded in 

L°°(0,T;H*(0,L)). 
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Figure 3: A sketch of mat. 



Moreover, for T small enough, we have 

< i? min < R + m(t, z) < iCax, WN eN,ze(0,L),te (o, T). (57) 



Proof. From Lemma 5.1 we have that < C, where C is independent of 
N. This implies 

\\VN(t)\\l>( ,L), \\9 Z VN(t)\\h(0,L), \\d 2 zz VN(t)\\l*(o,L) < C Vt G [0,T]. 
Therefore, 

ll^v|L°°(0,T;2£2(0,£)) < C- 

To show that the radius R + r/jy is uniformly bounded away from zero for T 
small enough, we first notice that the above inequality implies 

ll^-^llfl?(o,i)<2C, n=l,..., TV, NeN. 

Furthermore, we calculate 

n-l 



n-l 



WN-VO 1 1 £3(0,2,) 



< 



^at||l 2 (o,l) 



i+7 



% l|L 2 (0,L), 



i=0 



i=0 



5.1 



we have that E N 2 < C 



where we recall that rf N = tjq. From Lemma 
where C is independent of N. This combined with the above inequality 
implies 



\VN-Vo\\mo,L)<CnAt<CT, n = 1, . . . , N, N e N. 
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Now, we have uniform bounds for \\rf}^ — ^o||l 2 (o,l) an d ||*7jv — Vo\\h§(o,l)- 
Therefore, we can use the interpolation inequality for Sobolev spaces (see for 
example [1], Thm. 4.17, p. 79) to get 



VN-Vo\\m(o,L)<2CVT, n = 1, . . . , N, N G N. 



From Lemma 5.1 we see that C depends on T through the norms of the 
inlet and outlet data in such a way that C is an increasing function of T. 
Therefore by choosing T small, we can make ||^ — ^olltf^o,!,) arbitrary small 

for n = 1, , N, N G N. Because of the Sobolev embedding of H x (0, L) 

into C[0, L] we can also make \\r]^ — ?yo||c[o,L] arbitrary small. Since the initial 
data r/o is such, that R + r) (z) > (due to the conditions listed in ([8])), we 
see that for a T > small enough, there exist -R m i n , -R ma x > 0, such that 

< Rmin <R + 7] N (t, Z) < iCax, ViV G N, Z G (0, L) , t G (0, T). 

□ 

We will show in the end that our existence result holds not only locally 
in time, i.e., for small T > 0, but rather, it can be extended all the way until 
either T = oo, or until the lateral walls of the channel touch each other. 

From this Proposition we see that the L 2 -norm H/H^m) — / f 2 , and the 
weighted L 2 -norm H/H^a/™ = f(R + VN)f 2 ar e equivalent. More precisely, 
for every / G L 2 (Q), there exist constants C\,Ci > 0, which depend only on 
-Rmin, -Rmax, and not on / or N, such that 

c x [ (R+ VN )f 2 < WfWhm < C2 [ (R + VN)f 2 - (58) 

Jn Jq 

We will be using this property in the next section to prove strong convergence 
of approximate functions. 

Next we show that the sequences of approximate solutions for the velocity 
and its trace on the lateral boundary, are uniformly bounded. 

Proposition 6.2. The following statements hold: 

1- {vn)u&n is uniformly bounded in L°°(0, T; L 2 (0, L))flL 2 (0, T; Hq(0, L)). 

2. (u^-)neN is uniformly bounded in L°°(0, T; L 2 (0, L)). 

3. (uAr) n6N is uniformly bounded m L°°(0,T; L 2 {Q)) n L 2 (0,T; H\Q)) . 



30 



Proof. The uniform boundedness of (vjf) NeN, (v^) N<aNj an d the uniform bound- 
edness of (un)n £ fs in L°°(0,T; L 2 (f2)) follow directly from Statements 1 and 



2 of Lemma 5.1, and from the definition of (i>jv)neN, (vn)n£N and (un)ngn 



as step-functions in t so that 

f T 
Jo 



N-1 



n=0 



To show uniform boundedness of (u/v)jveN in L 2 (0, T; if we need to 



explore the boundedness of (Vujy)jveN- From Lemma 5.1 we only know that 



the symmetrized gradient is bounded in the following way: 

N 



£ 

n=l 



/ (R + V ^)\n%-\u n N )\ 2 At <C. (59) 
Jn 



We cannot immediately apply Korn's inequality since estimate (59) is given 
in terms of the transformed symmetrized gradient. Thus, there are some 
technical difficulties that need to be overcome due to the fact that our prob- 
lem is defined on a sequence of moving domains, and we would like to obtain 
a uniform bound for the gradient (Vu^r) NeN . To get around this difficulty we 
take the following approach. We first transform the problem back into the 
original domain fi^n-i on which is defined, and apply the Korn inequality 
in the usual way. However, since the Korn constant depends on the domain, 
we will need a result which provides a universal Korn constant, independent 
of the family of domains under consideration. Indeed, a result of this kind 
was obtained in [TTJ, HH], assuming certain domain regularity, which, as we 
show below, holds for our case due to the regularity and uniform boundedness 
of (^^T^TveN- Details are presented next. 

For each fixed JV 6 N, and for all n = 1, . . . , N, transform the function 
u^r back to the original domain which, at time step n, is determined by the 
location of the boundary r/ N at time step n — 1, i.e., by rf^f 1 : 

u(n, N) = u n N o A v n- h n = 1, . . . , JV, iV e N. 



By using formula (24) we get 



( J R + 7 7 -- 1 )|D^" 1 (^)| 2 = / |D(u(n,iV))| 2 = ||D(u(n, iV))||| 2(n n _ 
i Jn „_i v n 
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We can now apply Korn's inequality on f^n-i to get 
IVu^AOU^ <C( V ^ l )\\B(n(j,N))\\ 2 L2{n n _ l); n = 1, . . . , N, N E N, 



N 



where C{rf^~ ) is the Korn's constant associated with domain fi^n-i. Next, 
we transform everything back to Q by using the inverse mapping, and employ 
(58) to obtain: 

llV^u^lU^^^-^HDj" 1 ^)!!^^), n = l,...,N, NeN. 

-_n— 1 

Now, on the left hand-side we still have the transformed gradient Vjy and 
not V, and so we employ (25) to calculate the relationship between the two: 

Vu; = (V^uv) (vA,»->) , n — 1, . . . ,N, NeN. 

Since t]n are bounded in L°°(0, T; H 2 (0, L)), the gradient of the ALE mapping 
is bounded: 

HV^n-x \\ LOO(Q) < C, n = 1, . . . , N, N e N. 

Using this estimate, and by summing from n = 1, . . . , N, we obtain the 
following estimate for Vu^: 

N N 
n=l n=l J Q 

If we could show that 

C{r, n N l ) <K, n = l,...,N, N e N, 

we would have proved the uniform boundedness of in L 2 (0, T; H 1 ^)). 
The existence of a uniform Korn constant K follows from a result by Velcic 
p?] Lemma 1, Remark 6, summarized here in the following proposition: 

Proposition 6.3. [49J Let Q s C 1R 2 be a family of open, bounded sets with 
Lipschitz boundaries. Furthermore, let us assume that the sets Q s are such 
that fl s = F S (Q), where F s is a family of bi-Lipschitz mappings whose bi- 
Lipschitz constants of F s and F~ l are uniform in s, and such that the family 
{F s } is strongly compact in W 1,00 (Q 



,M 2 ) 
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Let u s G L 2 (Q S ,M. 2 ) be such that the symmetrized gradient D(u s ) = 
| (Vu s + (Vu s ) r ) is in L 2 (f2 s , IR 2 ). Then, there exists a constant K > 0, 
independent of s, such that 



| u s||vy 1 '2(n s ;K 2 ) < K 



udxidx^ 



(xi\i2 — X2Ui)dxidx^ 



\D(u) 



(60) 

We apply this result to our problem by recalling that Statement 1 in 
Lemma 5.1| implies 



m(o,L) <C, n = l,...,N, iVGN. 



Because of the compactness of the embedding H (0, L) CC W ,oo (0,L) and 



the definition of A j-i given in (22), we have 



||A.«-i||vri,ao < C, WA'^Wwi,*, <C,n = l,...,N, N e N. 
Furthermore, the set {A^n-i : n — 1, . . . , N, N G N} is relatively compact in 



W l '°°(Q). Thus, by Proposition 6.3 there exists a universal Korn constant 
K > such that 



N 



N 



El|Vu^||i 2(n) At<^^ / (l+ V n N - l )\B%-\u- N )\ 2 At, 

n=l n=l J Q 

which implies that the sequence (Vuat) n^n is is uniformly bounded in L 2 ((0, T) x 
Q), and so the sequence (u^n^ is is uniformly bounded in L 2 (0,T; if 1 (fi)). 

□ 

We remark that instead of Velcic's result, we could have also used the 
result by Chambolle et al. in [TT] . Lemma 6, pg. 377, which is somewhat 
less general, and in which the main ingredient of the proof is the fact that 
the fluid velocity is divergence free, and that the displacement of the domain 
boundary is only in the radial (vertical) direction. 

From the uniform boundedness of approximate sequences we can now 
conclude that for each approximate solution sequence there exists a subse- 
quence which, with a slight abuse of notation, we denote the same way as 
the original sequence, and which converges weakly, or weakly*, depending on 
the function space. More precisely, we have the following result. 
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Lemma 6.1. (Weak and weak* convergence results) There exist sub- 
sequences (^Ar)ATgN, (wnOatsn, (v%)n£N, and (un)n£N, and the functions r\ G 
L°°(0, T; H 2 (0, L)),ve L°°(0, T; L 2 (0, L))fU 2 (0, T; if 2 (0, L)), v* G L°°(0, T; L 2 (0, L)) ; 
and u G L°°(0,T; L 2 (n)) f] L 2 (0,T; H 1 ^)), such that 



7] N r] weakly* in L°°(0, T; # 2 (0, L)), 
v N v weakly in L 2 (0, T; # 2 (0, L)), 
v N v weakly* in L°°(0, T; L 2 (0, L)), 

v* N ->> weakly* in L°°(0,T;L 2 (0,L)), 

ujv -»> u weakly* in L°°(0, T; L 2 (fi)), 

u w ->> u weakly in L 2 (0, T; if 1 ^)). 



(61) 



Furthermore, 



(62) 



Proof. The only thing left to show is that v = v*. To show this, we mul- 
tiply the second statement in Lemma 5.1| by At, and notice again that 
\\ v N\\h«o,T)x(o,L)) = At En=i \\ v N\\h(o,Ly This implies \\v N -v* N \\ L 2 mT)x{0jL)) < 
Cy/At, and we have that in the limit, as At — > 0, v = v*. □ 



6.1 Strong convergence of approximate sequences 

To show that the limits obtained in the previous Lemma satisfy the weak form 
of Problem 3.1 we will need to show that our sequences converge strongly 
in the appropriate function spaces. To do that, we introduce the following 
notation which will be useful in the remainder of this manuscript: denote by 
Th the translation in time by h of a function / 

r h f(t,.) = f(t-h,.), ft, G R. (63) 



The strong convergence results will be achieved by using Corollary 6. 1| listed 
below, of the following compactness theorem [6]: 

Theorem 6.1. (Riesz-Prechet-Kolmogorov Theorem) Let f2 C M n be 

an open subset ofW 1 , and u C Q. Let J 7 be a bounded subset in L P (Q) with 
1 < p < oo. Assume that 

> 35 > 0, 5 < dist(w, <9fi), such that 

\\r h f - /\\lp( u ) <e,Vhe R n such that \h\ < 5 and V/ G 7. (64) 
Then F is relatively compact in L p (u). 
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Notice that u> C Q is introduced here so that the shifts f(x±h) would 
be well-defined. 

Corollary 6.1. (Corollary 4.37, p. 72 in [6j) Let Q C R n be an open subset, 
and let J 1 be a bounded subset in L P (Q), for 1 < p < oo. Assume that 



We > 0, Vw CC Q, 35 > 0, 5 < dist(w, <9fi), such that 
\T~hf - f\\Lp(u>) < e,\/h e W 1 such that \h\ < 5 and V/ G J 7 , 



and 



Ve > 3wCCfi such that 
Then, J 7 is relatively compact in L P (Q). 



LP 



(65) 



(66) 



The main ingredient in getting the "integral equicontinuity" estimate (65) 



is Lemma 5A Namely, if we multiply the third equality of Lemma 5.1| by 
At we get: 



\T~AtUN - u 7v|lL 2 ((0,T)xn) + ll r A^AT 



' ' U A r llL 2 ((0,T)x(0,L)) — CAt. 



(67) 



This is "almost" (65) except that in this estimate e depends on At (i.e., N), 



which is not sufficient to show equicontinuity (65). We need to show that 
estimate (65) holds for all the functions (vn)n£N, (u^Nm, independently 



of N G N. This is why we need to work a little harder to get the following 
compactness result. 

Theorem 6.2. Sequences (vn)n£N, (ujv)jv e N are relatively compact in 
L 2 (0,T;L 2 (0,L)) and L 2 (0, T; L 2 (£l)) respectively. 



Proof. The proof is based on Corollary 6.1 We start by showing that (65) 



holds. To do that, first notice that by Lemma 5^ both sequences are bounded 
(uniformly in iV G N) in the corresponding spaces. Furthermore, spatial 
derivatives (<9 2 t>Ar) n eN, (VuTvVeN are also bounded in L 2 (0, T; L 2 (0, L)) and 
L 2 (0, T; L 2 (fl)) respectively, which guarantees equicontinuity with respect to 
the spatial variables. Thus, to show that sequences (vn)n&, (utv) weN are 
"equicontinuous" in L 2 (0, T; L 2 (0, L)) and L 2 (0, T; L 2 (Q)), respectively, we 
only need to consider translations in time, r^. 



We proceed by proving "integral equicontinuity" (65 ) for sequence (vn) ngn- 
Relative compactness of (u/v)_/veN can be proved analogously. 

We recall, one 



Let e > and let C be the constant from Lemma 5.1 



more time, that C is independent on N, and thus of At. 
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Let u C f2 be an arbitrary compact subset of f2. Define 



5 := min{dist(u;, <90)/2, e/(2C)}. 

We will show that 

||T h Ujv - ^A r lli2( aJ; L2 (0iL)) < £, V|/i| < 5, independently of A/" e N. (68) 
Thus, for each JV 6 N, namely, for each At = T/N, we want to show that 



(68) holds, independently of N, for each h such that \h\ < 5. 

Let h be an arbitrary real number whose absolute value is less than S. We 



want to show that (68) holds for all At = T/N. This will be shown in two 



steps. First, we will show that (68) holds for the case when At > h (Case 
1), and then for the case when At < h (Case 2). 

A short remark is in order: For a given 5 > 0, we will have At < 5 for 
infinitely many N, and both cases will apply. For a finite number of functions 



(vn), we will, however, have that At > 5. For those functions (68) needs to 
be proved for all At such that \h\ < 5 < At, which falls into Case 1 bellow. 
Thus, Cases 1 and 2 cover all the possibilities. 



Case 1: At > h. We calculate the shift by h to obtain (see Figure |4j): 

N-l .4 At 



\\r h v N - v N \\ 2 L 2 {ui . L 2 {0}L)) < 2J / 

j=1 JjAt-h 



\ v n ~ IIl 2 (o,l) 



N-1 
3=1 

The last inequality follows from \h\ < 5 < e/(2C). 

Case 2: At < h. In this case we can write h = lAt + s for some I G N, 
< s < At. Similarly, as in the first case, we get (see Figure [B]): 



N-l-l 

Wv N ~ v N \\\ Hu] . L2 ^ L)) = ( /. 



^-'yl r (j+l)At-s 

\\ V N V N llL 2 (0,L) 

j=i ^' A * (69) 

+ / K-< +m lliw))- 

J(j+l)At-s 

Now we use the triangle inequality to bound each term under the two in- 
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f Case At>h 




CO 



Figure 4: Case 1: At > h. The graph of i>jy is shown in solid line, while the 
graph of the shifted function t^v^ is shown in the dashed line. The shaded 
area denotes the non-zero contributions to the norm Ht^at — ^tvII^- 



tegrals from above by Yli=i \\ v n 
terms together we obtain 



v 



j+i\\2 

n Hl 2 (o,l)- 



After combining the two 



N-l-l l+l 

- ^ \ V s II 2 ( 7 n.\ 

N ~ V N IIl 2 (0,LV l' U J 



\\r h v N - v N \\l2 {w . L 2 {0jL)) < At " 7 " 

j=l i=l 



Using Lemma 5.1 we get that the right hand-side of (70) is bounded by 



At(l + 1)C. Now, since h = I At + s we see that At < h/l, and so the right 
hand-side of (70) is bounded by l -j^hC. Since \h\ < 5 and from the form of 

5 we get 



l + l 



\r h v N - v N \\ L 2 {uJ . L 2 { ^ L)) < At(l + 1)C < -—hC < 



l + le 



I 2 



< 6. 



Thus, we have shown that (65) holds. 



To show that ((66]) holds, let e > 0. Define u = [e/(4C),T - e/(4C)]. 
We see that u is obviously compact in (0,T) and from the first inequality in 

Lemma 5.l| (boundedness of v ^r +2 , i — 1, 2 in L 2 (0, L)) we have 



(0,T)\w 
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f Case At<h=At+s 



t 



At 



to 



Figure 5: Case 2: At < h = At + s, < s < At. The graph of v N is 
shown in solid line, while the graph of the shifted function r^v at is shown in 
the dashed line. The shaded areas denote the non-zero contributions to the 
norm ||r fe t>7v — ^jvll^- The two colors represent the contributions to the first 
and second integral in ( 69 ) separately. 



By Corollary 6.1, the compactness result for (fjv)jveN follows. Similar 
arguments imply compactness of (ujv)jveN- □ 

To show compactness of (t]n)n€N we introduce a slightly different set 
of approximate functions of u, v, and 77. Namely, for each fixed At (or 
iV G N), define Ujv, rj N and vjy to be continuous, linear on each sub-interval 
[(n — l)At,nAt], and such that 

u N (nAt, .) = u N (nAt, .), v N (nAt, .) = v N (nAt, .), rj N (nAt, .) = n N (nAt, .), 

(71) 

where n — 0, . . . , N. See Figure [6} We now observe that 



d t f) N (t) 



7/ 



71+1 



7/'' 



At 



„n+l/2 _ „n 

At 



t e (nAt, (n + l)At), 



and so, since v* N was defined in (56) as a piece-wise constant function defined 
via v* N (t, •) = f n+ i, for t 6 (nAt, (n + l)At], we see that 

d t fj N = v* N a.e. on (0,T). (72) 



By using Lemma 



5.1 



(the boundedness of ~ : 



we get 



N)NeN 



is bounded in L°°(0, T; # 2 (0, L)) n W 1>oo (0, T; L 2 (0, L)). 
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We now use the following result on continuous embeddings: 

L°°{0,T;H%{0,L))r) W 1 ' oo (0,T;L 2 (0,L)) e — >■ C 0,1 - Q ([0,T]; # 2o (0, L)), (73) 

for < a < 1. This result follows from the standard Hilbert interpolation 
inequalities, see [H] . It was also used in [3T] to deal with a set of mollifying 
functions approximating a solution to a moving-boundary problem between 
a viscous fluid and an elastic plate. 

From (73) we see that (rj^NeN is also bounded (uniformly in N) in 
T];H 2a (0, L)). Now, from the continuous embedding of H 2a (0, L) 
into H 2a ~ e , and by applying the Arzela-Ascoli Theorem, we conclude that 
sequence (fj^Nm has a convergent subsequence, which we will again denote 
by (^at)at 6 n, such that 

fl N -> 77 in C([0, T]; # s (0, L)), < s < 2. 

Since sequences (r^TvgN and (r] N ) Ne ^ have the same limit, we have 77 = fj E 
C([0,T}; H s (0, L)), where r] is the weak* limit of (r] N ) N( z N , discussed in (61). 
Thus, we have 

fl N ^ r] in C([0,T];H s (0,L)), < s < 2. 

We can now prove the following Lemma: 

Lemma 6.2. rj N -»• r] in L°°(0,T; H s (0, L)), < s < 2. 

Proof. The proof follows from the continuity in time of rj, and from the fact 
that fj N ->■ 77 in C([0, T]; F s (0, L)), < s < 2. Namely, let e > 0. Then, 
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from the continuity of 77 in time we have that there exists a St > such that 



||77(ti) - 7)(t 2 )\\B-(Q,L) < 2? for ti,t 2 e [0,T], and \t 1 -t 2 \<8t. 

Furthermore, from the convergence fjj^ — > r\ in C([0, T]; H s (0, L)), < s < 2, 
we know that there exists an iV* G N such that 

\\Vn ~ v\\c{[o,T\;H s (o,L)) < |, VJV > iV*. 

Now, let N be any natural number such that iV > max{iV*, T/St}. Denote 
by At = T/N, and let t G [0, T\. Furthermore, let n G N be such that 
(n — l)Ai < t < nAt. Recall that fj^inAt) = t]N(nAt) = r)N(t) from the 
definition of t^jv and By using this, and by combining the two estimates 
above, we get 

\\r} N (t) - r}(t)\ I h*(o,l) = \\m(t) -r]{nAt) + i](nAt) - r](t)\\ H s (0>L) 

= \\r] N (nAt) -T](nAt) + 77 (nAt) - r)(t)\\ H s(o,L) 

< \\rj N (nAt) - r)(nAt)\\ + \\ V (nAt) - r)(t)\\ H s (0 , L) 

= \\fj N (nAt) - T](nAt)\\ H s {0:L) + \\rj(nAt) - r](t)\\ H s {0tL) < e. 

Here, the first term is bounded by e/2 due to the convergence rj^ — > i], while 
the second term is bounded by e/2 due to the continuity of r\. Since the 
obtained estimate is uniform in N and t, the statement of the Lemma is 
proved. □ 



We summarize the strong convergence results obtained in Theorem 6.2 



and Lemma 6.2 We have shown that there exist subsequences (un] 



NeN, 



(vn)n£N and (tj^n^n such that 



u N u in L 2 (0,T;L 2 (ft)), 

v N ->• v in L 2 (0,T;L 2 (0,L)), 

r At Ujv -> winL 2 (0,T;L 2 (ft)), (74) 

r At Vjv -> u in L 2 (0,T;L 2 (0,L)), 

t/jv -> « in L°°(0,T;if s (0,L)), < s < 2. 

Because of the uniqueness of derivatives, we also have v = d t r] in the sense of 
distributions. The statements about convergence of (rAtU^) n<=n an d (t^vn) tvsn 



follow directly from (67). 
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Furthermore, one can also show that subsequences (v N ) N and (un)n also 
converge to v and u, respectively. More precisely, 



u N u in L 2 (0,T;L 2 (fi)), 
v N -> v in L 2 (0,T;L 2 (0,L)) 



This statement follows directly from the following inequalities (see |47| . p. 
328) 



K - ^1112(0,^(0,^) < -IT l|r " ' ~ '' 



At A 

11 n+l „,n||2 



"TT / A\ v ~ v IIL2(0,L)> 



n=l 



At ^ 



N 

n+l ,,n||2 



,, ,i z << \ l, J- ,," z 



ra=l 



and Lemma 5A which provides uniform boundedness of the sums on the right 
hand-sides of the inequalities. 

We conclude this section by showing one last convergence result that will 
be used in the next section to prove that the limiting functions satisfy the 
weak formulation of the FSI problem. Namely, we want to show that 



The first statement is a direct consequence of Lemma [672] in which we proved 

3 
2 



that t]n — y Tj in L°°(0, T; H s (0, L)), < s < 2. This means that for s > § we 



immediately have 

■qN^rimL^T-C^L)). (77) 

To show convergence of the shifted displacements T^ t rj N to the same limiting 
function rj, we recall that 

fj N -> r/inC([0,r] ;j ff s [0,L]), < s < 2, 

and that {fj N ) Nm is uniformly bounded in C 0,1 ~ a ([0, T};H 2a (0, L)), 0<a< 
1. Uniform boundeness of (??Ar)/veN in C 0,1 ~ a ([0,T]; H 2a (0,L)) implies that 
there exists a constant C > 0, independent of N, such that 

\\fj N ((n - l)At) - rj N (nAt)\\ H 2 a{0tL) < C^t^. 

This means that for each e > 0, there exists an N\ > such that 

||7fcv((n - I) At) - ^(nAt)|| H2 « (0iL) < |,for all N > 
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Here, N% is chosen by recalling that At = T/N, and so the right hand-side 
implies that we want an N\ such that 

C (—} < - for all N > JVi. 
\N J 2 ~ 

Now, convergence fj^ — > ?y in CQO, T]; iP[0, L]), < s < 2, implies that for 
each £ > 0, there exists an N 2 > such that 

||^(nAt) - »7(*)l|jT-(o,i) < |, for all A^ > A^ 2 . 

We will use this to show that for each e > there exists an N* > max{A^i, iV 2 }, 
such that 

\\r A m(t) - v(t)\\H°(o,L) < e, for all N > N*. 

Let t G (0, T). Then there exists an n such that t G ((n — l)At,nAt]. We 
calculate 

\\TAtVN(t) -7?(*)||fP(0,L) = \\TAtVN{t) ~ VN^At) + fj N (nAt) - Tj(t) \\ H '(0fi) 

= \\f)N((n - I) At) - fj N (nAt) + fj N (nAt) - T)(t)\\ H *(o,L) 
< \\VN((n - l)At) - ^jv(nAt)||H«(o,L) + \\VN{nAt) - i](t)\\ H s {0tL) . 

The first term is less than e for all iV > N* by the uniform boundeness of 
(Vn) weN m C 0,1 ~ a ([0, T]; H 2a (0, L)), while the second term is less than £ for 
all N > N* by the convergence of fj N to 77 in C([0, T]; _£P[0, L]), < s < 2. 

Now, we notice that TAtVN — ( T AtVN)- We use the same argument as in 

Lemma 4.1. to show that sequences (ta^n) an d TAtVN both converge to the 
same limit r] in L°°(0,T; H s (0,L)), for s < 2. 



7 The limiting problem and weak solution 



Next we want to show that the limiting functions satisfy the weak form (47) 



of Problem |2.1| In this vein, one of the things that needs to be considered is 
what happens in the limit as A" — > oo, i.e., as At — > 0, of problem (53). Before 



we pass to the limit we must observe that, unfortunately, the velocity test 



functions in (53) depend of AH More precisely, they depend on because of 
the requirement that the transformed divergence-free condition • q = 
must be satisfied. This is a consequence of the fact that we mapped our 
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problem onto a fixed domain f2. Therefore we will need to take a special care 
in constructing the suitable velocity test function and passing to the limit in 



(53). 



7.1 Construction of the appropriate test functions 

We begin by recalling that the test functions (q, ip) for the limiting problem 



are defined by the space Q, given in (41 ), which depends on r\. Similarly, the 



test spaces for the approximate problems depend on N through the depen- 
dence on r] N . The fact that the velocity test function depend on N presents a 
technical difficulty when passing to the limit, as N — > oo. To get around this 
difficulty, we will restrict ourselves to a dense subset X of all test functions 
in Q , which is independent of rjN even for the approximate problems. The 
set X will consist of the test functions (q, ip) £ X — Xp x X$, such that the 
velocity components q of the test functions are smooth, independent of N, 
and V ■ q = 0. 

To construct the set Xp we follow ideas similar to those used in [TT] . We 
look for the functions q which can be written as an algebraic sum of the 
functions qo, which have compact support in U T in U T out U r&, plus a 
function q 1; which captures the behavior of the solution at the boundary 77. 
More precisely, let fl m i n and fl ma x denote the fluid domains associated with 
the radii R m in and R ma x, respectively. 

1. Definition of test functions (qo,0) on (0, T) x VL max : Consider all 
smooth functions q with compact support in fl v U T in U Y out U T b , and 
such that V • q = 0. Then we can extend q by to a divergence-free 
vector field on (0, T) x VL max . This defines qo- 

Notice that since rj^ converge uniformly to 7], there exists an N q > 
such that supp(qo) C ^r At ?? N , ViV > N g . Therefore, qo is well defined 
on infinitely many approximate domains n,. 4tW . 

2. Definition of test functions (qi,^) on (0, T) x VL max : Consider 
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V>eCi([0,T);# 2 (rg). Define 



qi 



A constant extension in the vertical 
direction of ipe T on Y v : qi := (0, ip(z)) T ; 
Notice divqx = 0. 



On Q m ax \ 



A divergence — free extension to Q r 
(see, e.g. [2E], p. 127). 



on fL 



From the construction it is clear that qi is also defined on fi TAtr)JV for 
each N, and so it can be mapped onto the reference domain Q by the 
transformation A TAtVN . 

For any test function (q, ip) G Q it is easy to see that the velocity component 
q can then be written as q = q — qi + qi, where q — qi can be approximated 
by divergence- free functions q , which have compact support in Q v U Y in U 
T out U T b . Therefore, one can easily see that functions (q, if}) = (q + qi,^) 
in X satisfy the following properties: 

• X is dense in the space Q of all test functions defined on the physical, 
moving domain VL V , defined by (41); furthermore, V ■ q = 0,Vq G Xp. 

• Similarly, for each q G Xp, define 

q = qoA rj . 

The set {(q, ^)|q = q o A^, q G Xp, ?p ^ ^s} is dense in the space Q v 
of all test functions defined on the fixed, reference domain Q, defined 



by (|46j). 

For each q G Xp, define 



q^v := q o A 



TAtVN- 



Functions q^r satisfy V rAtr?iV • qN = 0. These will serve as the test functions 
for approximate problems, defined on the domains determined by ta^at- The 
approximate test functions have the following uniform convergence proper- 
ties. 

Lemma 7.1. For every (q, tp) G X we have 

qiv - > q and VqAr — > Vq uniformly on [0, T] x Q. 
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Proof. By the Mean- Value Theorem we get: 



\qif(t,z,r) - q(t,z,r)\ = |q(*, z, (R + r At rj N )r) - q(t, z, (R + T])r)\ 

= \d r q(t,z,C)r\ \n(t, z) - rj N (t - At, z)\. 

The uniform convergence of q^ follows from the uniform convergence of 
since q are smooth. 

To show the uniform convergence of the gradients, one can use the chain 
rule to calculate 

d z q N (t, z,r) = d z q(t, z, (R+r At n N )r)+[d z r At n N (t, z)r] [d r q(t, z, (R + r At n N )r)] . 

The uniform convergence of d z qN{t, z,r) follows from the uniform conver- 
gence of Q z T A tT]N ■ Combined with the fist part of the proof we get <9 z qjv — > 
d z q uniformly on [0, T] x Q. The uniform convergence of d r q^ can be shown 
in a similar way. □ 
Before we can pass to the limit in the weak formulation of the approximate 
problems, there is one more useful observation that we need. Namely, notice 
that although q are smooth functions both in the spatial variables and in 
time, the functions q^ are discontinuous at nAt because T At rj N is a step 
function in time. As we shall see below, it will be useful to approximate each 
discontinuous function q^ in time by a piece- wise constant function, q^y, so 
that 

q N (t, .) = q(nAt-, .), te[(n- l)At, nAt), n = l,...,N, 

where qwinAt— ) is the limit from the left of q^ at nAt, n = 1, . . . , N. By 
combining Lemma 7_A with the argument in the proof of Lemma 62 , we get 

q^v — > q uniformly on [0, T] x Q. 



7.2 Passing to the limit 

To get to the weak formulation of the coupled problem, take the test functions 



(q N (t),ip(t)) (where q N = qoA 



T&tVN 1 



qeX F 



in equation (53) and integrate 
with respect to t from nAt to (n+ l)Ai. Furthermore, take ip(t) G Xg as the 
test functions in (51), and again integrate over the same time interval. Add 

,JV-1 



the two equations together, and take the sum from n — 0, 



to get 
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the time integrals over (0, T) as follows: 



Pf 



j j \r + T At 7] N ) (d t u N ■ q N + ^(r At u N - Wjv) • V TAtVN u N ■ q N 



~7:(r At u N - Wat) ■ V TAtm q N ■ u N 



P£ 
2 



T 

o Jn 



v N u N ■ q N 



+ 



o Jn 



(R + r A ^)2/iD^^(u N ) : D^( qjv ) + p s h 



T r L 



JO 



d t V N ll) 



with 



(a s (r) N ,il)) + a s (v N ,ip)) 

R([ P? n dt [ q z (t,0,r)dr- [ P Q N ut dt [ q z (t, L,r)dr) : 
Jo Jo Jo Jo 

V T ^. U ^ = 0, v N = ((u r ) N ) lr , 
u N (0, .) = u , 7/(0, .)n = Vo, v N (0, .) = v . 



(78) 
(79) 



Here Un and Wat are the piecewise linear functions defined in (71), r At is the 
shift in time by At to the left, defined in (63), V TAt,7JV is the transformed 
gradient via the ALE mapping A TAtVN , defined in (24), and v^, u;v, vn and 
t]n are defined in (56). 

Using the convergence results obtained for the approximate functions in 
Section |6j and the convergence results just obtained for the test functions 
qAr, we can pass to the limit directly in all the terms except in the term that 
contains dtUN- To deal with this term we notice that, since are smooth 
on sub-intervals (j'Ai, (j + I) At), we can use integration by parts on these 
sub-intervals to obtain: 

l-T r N-l Hi+l)At f. 

I (R + T~AtVN)d t u N ■ quv = / ( R + rf N )dtu N ■ 

Jo Jn J=0 JjAt Jn 

N-l r(j+l)At n 

= /A - I ( R + T AtVN)u N ■ d t q N 

j= o v JjAt Jn 
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+ / ( J R+77 J ' +1 -7f' +1 +^)u^ 1 -q JV (0' + l)At-)- / (R+v j W N -q N (jAt+)). 
Jn Jn ' 

(80) 

Here, we have denoted by qjv((j + l)At— ) and qjv(j'At+) the limits from 
the left and right, respectively, of q^v at the appropriate points. 

The integral involving dtqN can be simplified by recalling that q^v = 
q o A m , where 7]n are constant on each sub-interval (J At, (j + I) At). Thus, 
by the chain rule, we see that d t qN = d t q on (J At, (j + l)At). After summing 
over all j = 0, iV — 1 we obtain 

~ E / ( R + r AtW)ujv ■ d t q N = - (R + r A tVN)u N ■ d t q. 

• JjAt Jn Jo Jn 



To deal with the last two terms in (80) we calculate 

N-l 



E ( / (^ +1 -^ +1 +^)u; +1 -q^((j+l)At-)- / (R+rf N )u i N .q N (jAt+)) 
„_n y Jn Jn 7 



3=0 

N-l 



E / ((«+ ■ i*(0" + !) A *-) - - ^X 1 ■ <uv(0" + !) A *-)) 

„_n 

iV-1 

- / (i? + r/ )u • q(0) - V / (R + VnX ' <bfU*t+) 
Jn =1 Jn 



Now, we can write (r]^ +1 —rji) as f J+2 At, and rewrite the summation indexes 
in the first term to obtain that above expression is equal to 

= E / ( R + r l N ) u N ■ ^NijAt-) - f f v* N u N -q N 
• =1 Jn Jo Jn 



- f(R + Vo)u ■ q(0) - E / ( R + VnX ■ <lN(jAt+) 
Jn - =1 Jn 

Since the test functions have compact support in [0, T), the value of the first 
term at j = N is zero, and so we can combine the two sums to obtain 

N 

= E / ( R + VnX ■ MJAt-) - q N (jAt+)) 

3=1 Jn 
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- / (R + 7]o)u ■ q(0) - / / v* N u N ■ q N . 
Jn Jo Jn 

Now we know how to pass to the limit in all the terms expect the first one. 

We continue to rewrite the first expression by using the Mean Value Theorem 

to obtain: 

q N (jAt-, z, r)-q N (jAt+, z, r) = q(jAt, z, (R+rf N )r)-q(jAt, z, (R+rf^ l )r) = 

= d r q(jAt, z, ()r(rf N - r/^ +1 ) = -Atd r q(jAt, z, Qv^ 2 r. 
Therefore we have: 

N ~ 1 p i-T—At p 

V / {R+rf N )tf N (q{jAt-)-q{jAt+)) = - / (R+ VN )u N rT. At v* N d r q. 

• =1 Jn Jo Jn 

We can now pass to the limit in this last term to obtain: 

T-At r r T r 

(R + r] N )u N rT_ At v* N d r q (R + rj)\irdtr]d r q. 



Jo Jn 

Therefore, by noticing that <9 t q = dtq + rdtr]d r q we have finally obtained 

/ (R + r At r] N )d t u N ■ qjv -»> - / (R + T))u-d t q- / d t r)U ■ q 
o Jn Jo Jn Jo Jn 

- / (fl + 77o)uo-q(0), 
Jn 

where we recall that q = q o A v . 

Thus, we have shown that the limiting functions u and 77 satisfy the weak 
form of Problem 2.1 in the sense of Definition 4.2 , for all test functions (q, t/>), 
which are dense in the test space Q v . The following theorem holds: 

Theorem 7.1. (Main Theorem-The Viscoelastic Case) Let Qf, q s , fx, 

h, Ci, Di > 0, i = 1,2,3. Suppose that the initial data vq G L 2 (0,L), 
u G L 2 (Q Vo ), and i] e Hq(0,L) are such that (R + i] (z)) > 0, z e [0,L], 
Furthermore, let P in , P out G L / 2 OC (0, 00). 



Then there exist a T > and a weak solution of (u, rf) of Problem 3.1 



(or equivalently Problem 2.1) on (0,T) in the sense of Definition 1^.2 (or 
equivalently Definition J^.l), which satisfy the following energy estimate: 

E(t)+ / D(r)dr<E + C(\\P m \\ ), te[0,T], (81) 

Jo 
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where C depends only on the coefficients in the problem, E is the kinetic 
energy of initial data, and E(t) and D(t) are given by 

— ~^\\ u \\L 2 (n v (t)) + ~^~\\ a tm\L?{0,L) 

+ ^11^11^+^11^11 

D(t) = ^l|D(u)ll£ W t))) + D \\d t r]\\l 2 ^ L) + D 1 \\d t d z r ] \\l 2{0jL) + D 2 \\d t d 2 zV \\l 2{0tL) . 



Furthermore, one of the following is true: either 

1. T — oo, or 

2. lim min (R + n(z)) = 0. 

t^Tze[0,L] 



Proof. It only remains to prove the last assertion, which states that our result 
is either global in time, or, in case the walls of the cylinder touch each other, 
our existence result holds until the time of touching. This can be proved by 
using similar argument as in [TT] p. 397-398. For the sake of completeness 
we present the arguments here. 

Let (0, Ti), Ti > 0, be the interval on which we have constructed our solu- 
tion (u, rj), and let mi = mm( 0i Ti)x(o,L)(-R+ r ?)- From Lemma 6.1 we know that 



m l > 0. Furthermore, since 77 G W 1,DO (0, T; L 2 (0, L)) fl L°°[0, T; Hq(0, L)) 
and u G L°°(0,T;L 2 (n)), we can take T x such that t?(Ti) G H 2 (Q,L), 
d t r](Ti) G L 2 (0,L) and u(Xi) G L 2 (Q). We can now use the first part of 



Theorem 7.1 to prolong the solution (u, rj) to the interval (0,T 2 ), T 2 > T%. 
By iteration, we can continue the construction of our solution to the in- 
terval (0, T k ), k G N, where (T fc ) fcgN is an increasing sequence. We set 
m k = min (0 ,T fc )x(o,L)(# + v) > 0. 

Since > we can continue the construction further. Without loss 
of generality we could choose a T k+ i > T k so that m k+ i > From (81) 



and (73), by taking a = 1/2, we have that the displacement rj is Holder 



continuous in time, namely, 

IMIc .l/2(0,T fe+1 ;C[0,L]) - C(T k+ i). 

Therefore, the following estimate holds: 

R+r){T k+1 ,z) > R+ V (T k ,z)-C(T k+1 )(T k+1 -T k y/ 2 > nik-CiT^T^-T^ 1 / 2 . 
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For a T k+ i chosen so that m k+ i > ^ this estimate implies 



m 



T k+1 -T k > A ^ f J ken. (82) 



Now, let us take T* = sup km T k and set m* = min( 0j T*)x(o,L)(-R + v)- 
Obviously, m k > m*, k G N. There are two possibilities. Either m* = 0, 
or m* > 0. If m* = 0, this means that lim min (R + n(z)) = 0, and the 

t-+T z e[0,L} 

second statement in the theorem is proved. If m* > 0, we need to show that 
T* = oo. To do that, notice that (81) gives the form of the constant C(T) 
which is a non-decreasing function of T. Therefore, we have C(T k ) < C(T*), 
\/k G N. Using this observation and that fact that m k > m*, k G N, estimate 
(82) implies 

(m*) 2 

n + i-n>^^y 2 , v/cgn. 

Since this holds for all k G N, we have that T* = oo. □ 



8 The purely elastic case 

We consider problem with the structural viscoelasticity constants 

D = D\ = D 2 = 0. Thus, we study a FSI problem between the flow of a 
viscous, incompressible fluid, and the motion of a linearly elastic Koiter shell. 



We show here that an analogue of Theorem 7.1| holds, with a proof which is 
a small modification of the proof presented for the viscoelastic case. 

We begin by defining weak solutions. First notice that the following 
formal energy inequality holds, which is an analogue of the energy inequality 

j t E(t) + D F (t) < C(P m (t),Pout(t)). 

Here Dp{t) = //||D(u)||| 2 m includes fluid dissipation, while E(t) ac- 
counts for the kinetic energy of the fluid and of the structure, and the elastic 
energy of the Koiter shell. Thus, E(t) is the same as in the viscoelastic case, 



and is given by (15). 

To define a weak solution we introduce the solution spaces. The only 
difference with respect to the viscoelastic case is in the solution space for the 
structure, which has lower regularity in time. Namely, we set 

W SE (0, T) = W 1 ' 00 ^, T; L 2 (0, L)) n L°°(0, T; V 5 ), (83) 
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to replace the space Ws, defined in (39), where Vs = Hq(0, L) as before. The 
velocity space is the same as before, and is given by Wf, defined in (38). 
The solution space for the coupled FSI problem is now given by 

W E (0,T) = {(11,77) G Wf(0,T)xW se (0,T) : u(t, z, R+r]{t, z)) = d t r}{t, z)e r }, 

(84) 

which replaces the space (40) for the viscoelastic case. 

For the problem defined on a fixed, reference domain, the solution space 

W E (0, T) = {(u, n) G W£(0, T) x W SE (0, T) : u(t, z, 1) = d tV (t, z)e r } (85) 



replaces W v defined in (45). 

With these modifications, we have the following analogue definition of 
weak solutions: 

Definition 8.1. We say that (u, 77) G Wb(0, T) is a weak solution of Prob- 
lem \2l\ with D = D 1 = D 2 = 0, if for every (q,-0) G C^([0, T);Vp x V SE ) 
such that q(i, x, R + rj(t, x)) = i/j(t, x)e r! the following equality holds: 

Pf(~ [ [ u-«9 t q+ f 6(t,u,u,q))+2^ / / D(u):D(q) 

Jo Jn n (t) Jo Jo Jn„(t) 



Pi 
2 



T pL 



JO 



(d t r)) 2 ip - p s h 



T pL 



JO 



d t 7)d t i/) + / 05(77, i>) 



= [ (F(t), q) rm/out +p f f Uo • q(0) + p s h [ v ^{0). 
Jo Jn vo Jo 

(86) 

For the problem defined on a fixed, reference domain, the definition of 
a weak FSI solution for the purely elastic structure case is given by the 
following. 

Definition 8.2. We say that (11,77) ^#(0, T) is a weak solution of Problem 



3.1 with D = Di = D 2 = 0, if for every (q, if)) G C*([0,T); V F x V SE ) such 
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that q(t, z, 1) = if)(t, z)e r , the following equality holds: 

p f (- f [ (R + r ] )u-d t q+ A"( U ,u,q)) 
Jo Jn Jo 

+2p f [ (R + r))W(u) :D"(q)-^ f [ (d t r))u ■ q 
Jo Jn 1 Jo Jn 



-Psh [ [ d t r]d t ij+ [ a 3 (r],i>) (87) 
Jo Jo Jo 



= R (P in (t) / (q z )\ z=0 - Pout(t) / (q z )\z=L) 
Jo Jo Jo 

+p f u ■ q(O) + p s h / Wo^(O). 

J!J, JO 

Theorem 8.1. (Main Theorem - The Elastic Case) Let Qf, g s , p, h, 

Ci > 0, Di = 0, i = 1,2,3. Suppose that the initial data Vq G L 2 (0,L), 
u G L 2 (Q V0 ), and i] G Hq(0,L) are such that (R + i] (z)) > 0, z G [0,L]. 
Furthermore, let Pi n , P out G Lf oc (0, oo). 

Then, there exist a T > and a weak solution (11,77) on (Oj^") of Problem 



3.1 (or equivalently Problem 2.1) with D = D\ = D 2 = 0, in the sense of 



Definition 8.2 (or, equivalently, Definition 8.1), which satisfies the following 



energy estimate: 



E(t)+ [ L> F (r)dT<£; + C(||P, ri ||| 2(0!t) + ||P out ||i 2(0it) ), tG[0,T], (88) 
Jo 

where C depends only on the coefficients in the problem, Eq is the kinetic 
energy of initial data, and E(t) and Dp(t) are given by 



E(t) 



P±\ 

I 



u 



L 2 (n v (t)) 



, P sh \\A 112 

+ -Y\\dtV 



L 2 (0,L) 



+ ACo\\v\\l>(p,L) + Ci\\dzV\\h(o,L) + c 2\\d 2 z v\\h(o, L ))^ 



D F (t) = p\\B(u)\\ 2 L2 



(n„(t)))- 



Furthermore, one of the following is true: either 
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1. T = oo, or 



2. lim min (R + riiz)) = 0. 
t^Tze[o,L] y 



Proof. The main steps of the proof are analogous to the proof of Theorem 



7.1 We summarize the main steps and emphasize the main points where the 
proofs differ. 

The problem is again split into the structure elastodynamics problem, 
i.e., problem Al, and the fluid problem with the Robin-type boundary con- 
dition, i.e., problem A2. Problem Al is the same as in the viscoelastic case. 
Problem A2 is different since only the structure inertia (the structure accel- 
eration) is now coupled to the fluid motion via the boundary condition at 
the moving boundary, since Di = 0,i = 0,1,2. 

The semi-discretization of Problem Al is the same as before, including 
the analysis part in which the existence of a weak solution is proved. The 
semi-discretization of Problem A2 differs by the fact that now the bi-linear 
functional associated with structural viscosity a' s = 0. This influences the 
solution spaces: a weak solution to Problem A2 now belongs to: 

(u n+1 X +1 ) G Vjf x L 2 (0,L), (instead of Vf x # 2 (0,L)). 

The proof of the existence of a weak solution follows the same arguments 
based on the Lax-Milgram Lemma as before. 



Furthermore, an analogous energy estimate to (|55|) holds, with the 

n+ 
JV 



semi-discrete viscous dissipation -D^ +1 replaced by 



(D F ) n N +1 = Atfi [ (R + V n )\D*> n (u" N +1 )\ 2 . 
Jn 

Thus, up to this point, as before, we have a time-marching splitting 
scheme which defines an approximate solution to our FSI problem, in which 
the structure is modeled as a linearly elastic Koiter shell. Furthermore, for 
each At, the approximate FSI solution satisfies a discrete version of the en- 
ergy estimate for the continuous problem. 

The uni form energy bounds for the approximate solutions follow from 
Lemma [HTT] where (Dp) 71 ^ 1 replaces LP* 1 . The same method of proof pro- 
vides the statements 1-4 of Lemma 
the purely elastic case. 



5.1 



with (D F )^ +1 replaced by D 7 ^ 1 in 
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Therefore, we can proceed to get a weak convergence result, i.e., an 
analogue of Lemma |6.1| Of course, we will not have the weak convergence 
v n — v weakly in L 2 (0, T; Hq(0, L)), however, this convergence property was 
only used in passing to the limit in the bi-linear form a' s , which no longer 
exists in the purely elastic case, and is therefore not needed in the proof. 
Thus, a summary of the weak and weak* convergence results for the purely 
elastic case is given by the following: 

r] N ->> n weakly* in L°°(0, T; # 2 (0, L)), 

v N -± v weakly* in L°°(0, T; L 2 (0, L)), 

v* N ->> v* weakly* in L°°(0, T; L 2 (0, L)), (89) 

u N -»> u weakly* in L°°(0, T; L 2 (fi)), 

u N ->> u weakly in L 2 (0,T; iZ 1 ^)). 

The next step is to obtain a compactness result which will be used 
to show that the limiting functions satisfy the weak formulation of the FSI 



problem, i.e., we would like to obtain an analogue of Theorem 6.2 This 



is the place where the viscoelastic and purely elastic case differ the most. 



Namely, in the proof of Theorem 6.2 we used the fact that (<9 2 t>Ar) ne N is 



bounded in L 2 (0, T; L 2 (0, L)) (which followed from the uniform bound of the 
structural viscous dissipation terms). This allowed us to have equicontinu- 
ity with respect to the spatial variable of the sequence (fjv)neN m L 2 (0,L), 
and use Corollary |6.1| to obtain the corresponding compactness result in 
L 2 (0, T; L 2 (0, L)) by considering only translations in time to show an integral- 
type "equicontinuity" result in time for the sequence (wAr)neN- Since the vis- 
coelastic structural terms are no longer present, the regularization of due 
to structural viscosity can no longer be used. 

To get around this difficulty, we first recall that vn = (ujv)[r • e r . This 
follows from the definition of approximate functions ( }56| ) in which vn is 
defined as a solution of Step A2 for a given iV G N (see the comment below 



(56)). Using the Trace Theorem we then have 

ll V w|U 2 (0,T;# 1 /2(o,.L)) < C1|WJV|U 2 (0,T;/P(Q)), N G N, 

where C does not depend on N (it depends only on Q). Therefore, we have 
that (fAr)neN is uniformly bounded in L 2 (0, T; iZ" 1//2 (0, L)). Now, because of 
the compactness of the embedding if 1//2 (0, L) CC L 2 (0, L), the compactness 
in L 2 (0,L) of the sequence (fAr) ne N with respect to the spatial variable, fol- 
lows immediately. Therefore, similarly as before, we only need to consider 
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translations in time to prove integral "equicontinuity" in time, and use Corol- 
lary 6.1 to obtain that (fjv)neN is relatively compact in L 2 (0, T; L 2 (0, L)). 



This part of the proof is identical to that of Theorem |6.2 

Showing relative compactness of (u A r) neN is the same as in the viscoelastic 
case. 

Therefore, at this point we have the following strong convergence re- 
sults: 

U7V ->■ u in L 2 (0, T; L 2 (Q)), 

v N vmL 2 (0,T;L 2 (0,L)), 
r At u N -> umL 2 (0,T;L 2 (Q)), {JU> 
r At v N ->• u in L 2 (0,T;L 2 (0,L)). 

To show strong convergence of t]n we used only the information from the 
elastodynamics of the structure (Step Al), and the uniform bounds provided 



by Lemma 5.1 involving the kinetic energy of the fluid, the kinetic energy 



of the structure, and the structural elastic energy. See (72) and the proof 
below. All this is the same as in the purely elastic case, and so the same 
strong convergence result holds for rj^ as before: 

y Tj in L°°(0, T; H"(0, L)), < s < 2. (91) 

The rest of the proof is the same, since the construction of the appro- 
priate test functions only uses convergence properties of tjn, which are the 
same as in the viscoelastic case. 

Passing to the limit is performed in an analogous way as before. 

□ 

Therefore, slight modifications of the proof in the viscoelastic case, sum- 
marized above, lead to the proof of the existence of a weak solution to a FSI 
problem between the flow of a viscous, incompressible fluid, and a linearly 
elastic Koiter shell model, as stated in Theorem |8.1| This shows, among 
other things, the robustness of the method, presented in this manuscript. 



9 Appendix: The linearly viscoelastic Koiter 
shell model 

In this section we summarize the derivation of the linearly viscoelastic cylin- 
drical Koiter shell model, presented in equation Q, and list the expressions 
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for the coefficients Cj, Di, i = 0,1, 2, in terms of Youngs modulus of elasticity 
and Poisson ratio. More details of the derivation can be found in [TJ |9]. 

Consider a clamped cylindrical shell with the reference radius of the mid- 
dle surface equal to R > 0, with the shell thickness h > 0, and cylinder length 
L. The basic assumptions under which the Koiter shell model holds are: (1) 
the shell is thin (h/R << 1); (2) the strains are small, although large deflec- 
tions are admitted, and the strain energy per unit volume of the undeformed 
body is represented by a quadratic function of strain for an isotropic solid 
(Hooke's law); and (3) the state of stress is approximately plane. 

Denote by = (£ 2 (z), (, r ( z )) the displacement of the middle surface 
at z, where £ x (z) and £ r (z) denote the longitudinal and radial component 
of displacement respectively. Here, the axial symmetry of the problem has 
already been taken into account assuming that the displacement in the 9- 
direction is zero, and that nothing in the problem depends on 6. The change 
of metric and the change of curvature tensors for a cylindrical shell are given 
respectively by 



7«) 



t' z 

R£ r 



fr 



Here ' denotes the derivative with respect to the longitudinal variable z. 
Introduce the following functional space 

V c = Hl(0,L) x H 2 (0,L) = (fe,£ r ) G H\0,L) x H 2 (0,L) : 
6(0) = UL) = 6(0) = ML) = 0,£(0) = C(L) = 0} . 

Then a weak formulation of a linearly elastic cylindrical Koiter shell is given 
by the following: find r] = (r] z ,r] r ) G V c such that 



h f L h 3 
A 1 (rj)-^)Rdz + 



2 7 ' w/ 24 Jo 



Aq{t)) ■ Q{£)Rdz 



f-ZRdz, deV c , 



where ■ denotes the scalar product A ■ B := Tr (AB T ) , A,B G Af 2 (R) = M 4 . 
Here / is the surface density of the force applied to the shell, and A is the 
elasticity tensor given by [Hj: 



A, 



4A/i 
\ + 2^ 
1 
R 2 



A c ■ E) A c + 4 / uA c EA c , E G Sym i 



with 



A c 



1 
o Wl 
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where A and // are the Lame constants. Using the following relationships 
between the Lame constants and the Young's modulus of elasticity E and 
Poisson ratio a 

2/iA „ A + n E 2/xA H/i 1 A E 

\-2fi = 4/x = -, — — =4/i -— — = - -a, 



A + 2/i ^ ^A + 2/i 1-(T 2 ' A + 2,u A A + 2/i2A + /i 1 - a 
the elasticity tensor ^4 reads 

AE = -^^-(A c • E)A C + -^^A C EA C , E e Sym (M 2 ). 
1 — cH 1 + a 

Integration by parts in the weak formulation gives rise to the following static 
equilibrium equations in differential form: 




The Linearly Elastic Cylindrical Koiter Shell Model 

The terms multiplying h/2 account for the stored energy density due to 
stretching (membrane effects) and the terms multiplying h 3 /12 account for 
the stored energy density due to bending (flexural shell effects). 

To include the viscoelastic effects, we assume that rj is also a function 
of time. Viscoelasticity will be modeled by the Kelvin- Voigt model in which 
the total stress is linearly proportional to strain and to the time-derivative 
of strain. For this purpose we define the viscosity tensor B by: 

BE = ^^(A c • E)A C + -^A C EA C , E e Sym(M 2 ). 
1 — o v 1 + cr v 

Here E v and a v correspond to the viscous counterparts of the Young's mod- 
ulus E and Poisson's ratio a. Then, for a linearly viscoelastic Koiter shell 
model we define the internal (stretching) force 

N := ^7(77) + ^7(77), 



and bending moment 
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The weak formulation of the dynamic equilibrium problem for a linearly 
viscoelastic Koiter shell is then given by the following: for each t > find 
r](-,t) e V c such that V£ 6 V c 



+p s h J ■ iRdz = J f ■ £i?Gfe, 

where p s denotes the volume shell density. 

After integration by parts, the weak formulation above implies the fol- 
lowing dynamic equilibrium equations in differential form: 



d 2 r] z ~ drj r 
Psh dtl ° 2 dz 


3 dz 2 


jj d 2 rj r 
dtdz 


3 dtdz 2 


f Z1 


d TJr ~ 


r, d % 

dz 2 

1 dtdz 2 


+ 0^ + 

dz 


° A dz\ 






+ D Vz 
dtdz 


+ DA dtdz± 


fri 



The Viscoelastic Cylindrical Koiter Shell Model 



where 



~ hE h 2 ~ _ h 3 Ea - h Ea 



R 2 {l-a 2 y 12R 2/ ' 1 6R 2 (l-a 2 )' z Rl-a 2 
* _JiE_ ~ _ h 3 E 

u R 2 v l2R? h 6„ R 2 R 

h 3 

Ds = hC v , D4 = j^Cv, 



(92) 



and 



E E a 

- = ~, 5"? E) v .= - -. 
1 — (J 1 — <7 



Mathematical justification of the Koiter shell model can be found in [15| 

SUES ESI HQ]. 

In our problem we assume that longitudinal displacement is negligible, 
i.e., r] z = 0, as is commonly used in modeling blood flow through human 
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arteries. Thus, the resulting dynamics equilibrium equations in differential 
form are given by 



Psh-^ + C 0Vr ~C 1 -^ + C 2 ^ + D —- Dl — + D 2 — = f r 



where 



hE h 2 h 3 Ea h 3 E 



R 2 (l-a 2 ) y YlR? n 1 QR 2 {l-a 2 Y z 121-a 2 ' 
h h 2 h 3 D h 3 

This is exactly Q. 
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